General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 




I I5S5S5""P£»hh 

■■■■■■■■■■■■■■■■ 
■■■■■■■■■■■■■■■■ 




iNASA-Ca-17l20d) Diif I NIT II N Cf OfcCUND TEST 

FOU VEtilFICAl'lCh Li L AIM t SPACE SlliUCTUR E 

CONTROL fxLal lepoit, jQ Stp. 1980 - 01 

JUi. 1*d<* (control Oyraujict Cc.) 157 p 

iiC A08/.1F AO 1 CSC 1 2 GJ/18 


iiSiiiSSSSMiii 





DEFINITION OF GROUND TEST 
FOR VERIFICATION OF LARGE 
SPACE STRUCTURE CONTROL 


FINAL REPORT 
NOVEMBER 1984 


Sponsored By: 

George C. Marshall Space Flight Center 
Marshall Space Flight Center, Alabama 35812 

Under: 

Contract NO. NAS8-34700 


Contributors: 

Dr. George B. Doane III 
Dr. John R. Glaese 
Danny Tollison 
Thomas Howsman 
Sally Curtis 
Betsy Banks 

Edited By: 

Sally Curtis 


Prepared By: 


Control Dynamics Company 
555 Sparkman Drive, Suite 1414 
Huntsville, AL 35805 


ABSTRACT 


This paper discusses the work performed under contract No. NAS8-34700. 
The body of the paper deals mainly with descriptions of the tasks undertaken 
in a chronological order, with references made to detailed theories contained 
in the monthly reports. Those reports referenced are located in their corre- 
sponding appendices so that this document may stand alone. 

The main tasks discussed in this report include control theory and 
design, dynamic system modelling, and simulations of test scenarios. Other 
miscellaneous projects also in this paper include two papers presented and 
co-authored by Marshall Space Flight Center and Control Dynamics, and super- 
vision of hardware construction. 
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INTRODUCTION 


This contract, throughout the years of its existence, has been associ- 
ated with several different and varying tasks. These tasks have ranged from 
evaluating test articles and testing locations to control system design to 
system modeling to test simulations. These tasks and others all maintain the 
goal of supporting Marshall Space Flight Center (MFSC) in its effort to 
achieve a successful ground test experiment of a large space structure (LSS). 

This final report deals with the time period ranging from October 1982 
until the contract's official end in July 1984. A report dealing with the 
first year of the contract has previously been delivered in the form of an 
oral presentation to MSFC and can be found in Appendix A. 
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LEVEL of EFFORT 

The first project in the reporting time frame was the generation of a 
planar model of the ground test experiment structure. This was done by 
adapting an in-house modal analysis code (Gimbalflex) to allow placement of 
sensors and actuators at certain locations on the structure. At this point 
certain system characteristics had to be assumed as they were not as well 
known as they are currently. Figure 1 depicts the planar model and Attach- 
ment 1 from the 13th Monthly Report (MR) (found in Appendix B) delves more 
carefully into the theory behind this model. 

The next step in the development of this model was to adapt it for use 
with the control pole placement algorithm. This amounted to removing the 
uncontrollable rigid body mode from the model. This uncontrollable mode was 
due to the pointing of the gimbal . The method to eliminate this mode is 
discussed in Attachment 1 from MR #14 (Appendix C). 

Control system design for the planar model was also initiated during 
this time. A digital controller was being planned utilizing the Control Pole 
Placement technique. One assumpticn made in dealing with this design was the 
existence of rate of state feedback. 

Once this controller and other control algorithms were developed, it 
became necessary to conduct speed tests on the HP9845C to determine if the 
20ms sample period requirement could be met. It was discovered that it could 
not be met with the given equipment. A study was then made of additional 
hardware/ software necessary to meet the speed requirement. The results of 
this study yielded several alternatives varying in cost and can be found in 
the 16th MR (Appendix D). 
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EARTH FIXED (LABORATORY FRAME) 


FIGURE 1. Planar Version of the LSS-GTV Structure. 
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Next, an analysis was made concerning the support of the tip instrument 
package mount. The following recommendations were made: 

1. The tip support mechanism should be capable of exerting a constant, 
but adjustable force on the tip of the structure. The magnitude of 
this force should be variable over the range from a value corre- 
sponding to all the weight of the equ pment above the roll gimbal 
down to as near to zero as practical (most likely around 5 pounds). 

2. The mount used to attach the tip instrument package to the mast 
should incorporate the possibility of attaching weights to it. 

These weights could be used not only to balance the instrument pack- 
age with respect to the mast but also to cause known and adjustable 
degrees of unbalance. This unbalance could be used to couple iner- 
tially the different modes of oscillation; i.e. in plane, out of 
plane, and torsional . 

Figure 2 depicts the system design at the time of the above recommendations. 

In order to design an accurate control system, a good dynamic model is 
necessary. At this stage it was deemed necessary to construct a better model 
of the ASTROMAST than had previously existed. The beam was modeled as 273 
point masses (three for each of the 91 sections), each mass having three 
translational degrees of freedom with the exception of the three masses 
representing the tip section. They were restrained in one degree of freedom 
by being attached to the tip mass instrument package. This yielded a system 
of 816 second order equations of the form: 


Mx + Kx = 0, 
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where M Is a diagonal matrix of order 816 and K is a symmetric, 14 banded 
matrix, also of order 816. 

It can be easily seen that the problem of determining the modal fre- 
quencies of the system can be reduced to that of determining the values of 
w (frequency) such that 

det[ K - M w 2 ] - 0. 

This approach was straightforward; however, computational problems 
remained in finding the determinant of an 816th order matrix. Several 
methods were attempted which took into account the sparseness, bandedness, 
and symmetry of the problem. While these methods did reduce the problem 
somewhat, the problem was still extensive. 

At this point a less complex model was developed. This method Involved 
dividing the ASTROMAST into several lumped masses joined by massless rods. 
This Is a standard practice In dealings with beam structures as developed by 
N.O. tyklestad and M.A. Prohl . The number of lumped masses was varied, until 
the value of 91 was eventually decided upon. Each section was described by 
a point and field matrix relating deflections, slopes, moments, shears, 
angles, and torsions at each lumped mass. The matrices were then consecu- 
tively multiplied to relate base and tip conditions. Because of the symmetry 
of the ASTROMAST, the torsional and translational modes could be uncoupled to 
obtain these frequencies separately. The theory behind this method can be 
found in the 21st MR (Appendix E). The values obtained using this method are 
listed in Table 1. Independent checks for both torsional and translational 
frequencies were calculated and they compared favorably. (Torsional check 
may also be found in Appendix E.) 
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TABLE 1 

MODEL FREQUENCIES 


Bending Frequencies 

Torsional Frequencies 

0.62 Hz 

6.88 Hz 

3.92 Hz 


11.18 Hz 



Once this model of the mast was built, It could then be incorporated 
into a modal synthesis technique to build up a model of the entire test 
system. The results of this model are located in Appendix F. 

Using the above results, Dr. Henry Waites of MSFC and Dr. George Doane 
and Danny Toll Ison of Control Dynamics wrote and presented a paper at the 
AIAA Guidance and Control Conference at Gatllnourg, Tennessee during August 
of 1983. This paper was entitled: "Definition of Ground Test for Large 

Space Structure Control Verification". This paper described the experiment 
Itself, the equipment to be used, the above modal results and a comparison of 
these results with the experimental results. The comparison results are shown 
below in Table 2. A copy of this paper is located in ,‘ppendix G. 

At this point the work shifted back to control related tasks. A review 
was made of Dr. Henry Waites' technical papers on controlling large space 
structures. Three techniques were studied: 

1. observer techniques, 

2. closed-loop pole placement techniques, and 

3. a disturbance Isolation technique. 
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Experimental (Hz) 

Analytical (Hz) 

X Difference 

Bending* 

0 . 56 

0.62 

11 


3.4 

3.9 

lb 


— 

11 

— 

- 




Torsion 


6.9 



21 

21 

0.0 


29 

34 

19 

=2 


* Two modes at each frequency 


* 

Table 2. Summary cf Modal Frequencies from Model and Test Results for i,;.e 
Cantilivered ASTROMAST Beam. 
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Of these three methods, the Pole Placement Method is expected to be the most 
practical and useful. For a more thorough examination of these methods, see 
the 23rd MR (Appendix H). 

Another important aspect of this contract was simulation work. A simu- 
lation existed written in APL on the NASA computer, but one needed to exist 
at Control Dynamics in FORTRAN. This simulation would eventually be used to 
run test scenarios of the actual test set-up. The simulation contains trans- 
formation matrices from the tip and base instrument packages to the labora- 
tory reference frame, the strapdown algorithm, torque controllers. Dr. Waites' 
controllers, and the results from modal analysis techniques. 

The first step in this process was to convert the APL simulation to 
FORTRAN and resolve any discrepancies which might appear between the two ver- 
sions. Next, various details were added to achieve a more realistic test 
scenario; uncertainties were added to the gyros and accelerometers, and accel- 
erometer readings were added at the tip and base to allow for corrections due 
to constant accelerations due only to gravity. 

Dr. Waites of MSFC and Danny Tollison of Control Dynamics wrote and 
presented another paper entitled: "Ground Test Experiment for Large Space 
Structures" at the 16th Southeastern Symposium on System Theory at Mississippi 
State University. For more information on this paper see the 28th MR 
(Appendix I). A major modification initialed by MSFC, of the test structure 
appeared in this paper: the entire structure had been flipped up-side down. 
This can be seen in Figure 3. Several reasons existed which brought about 
this change; a main reason was the complexity and cost of the tip support 
system. The revised structure reduced this to a more manageable problem. 
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Again, the dynamics came into question. A better model was again needed 
not only for the control system, but also to be modified to obtain LSS 
characteristics (densely packed frequencies between .1 and .2 Hz). For this 
model the test structure was approximated as a continuous beam with several 
equivalent finite elements. 

At first, a model was made consisting only of the ASTROMAST in a canti- 
levered state. The first few frequencies (bending-torsional) obtained with 
this method agreed to within 5* of the experimental values shown in Table 3. 

Next, the remaining elements of the test structure were added to obtain 
a total model: shake table, AGS, roll gimbal , ASTROMAST and tip package. 

This model yielded five rigid body modes including pendulum modes at .15 Hz. 
Table 4 contains the next five modal frequencies. 

As can be seen from the tables, the frequencies are not densely packed. 
A couple of methods to try and achieve this characteristic were brought under 
consideration: 

1. Add weighted flexible stringers alongside the ASTROMAST. 

2. Add a flexible cruciform shape onto the bottom of the tip mass. 
Neither of these methods were implemented under this contract. However, the 
second method was used in the follow-on contract. 

Some remaining tasks to be dealt with on this contract included modify- 
ing the Work Breakdown Structure. As the project advanced and the testing 
phase approached, it was seen that changes to the WBS were necessary. These 
changes were initiated during this contract, but was not concluded until the 
next contract came into effect. 

Two hardware tasks existed on this contract. Control Dynamics was to 
oversee these tasks. The first task involved the construction of an air 
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bearing for the roll gimbal . This was deemed necessary when the structure 
was turned up-side-down, as the shaft in the roll gimbal was no longer 
restrained and could dislodge itself. The air bearing would hold it in 
place. The second task involved constructing a new tip instrument package 
bracket out of aluminum, as the previous bracket was made from steel and was 
too massive for its purpose. 


TABLE 3. CANTILEVERED ASTROMAST FREQUENCIES 


MODE 

MODEL 

VALUES (Hz) 

EXPERIMENTAL 

VALUES 

% DIFFERENCE 

1st BENDING 

0.543 

0.555 

-2.16 

2nd BENDING 

3.411 

3.430 

-.55 

1st TORSION 

4.016 

4.015 

.02 

3rd BENDING 

9.65 

8.83 

9.3 

2nd TORSION 

12.83 

12.23 

4.37 


TABLE 4. SYSTEM FREQUENCIES 


MODE # 

6 

7 

8 
9 

10 


FREQUENCY (Hz) 

1.12 

1.19 

1.34 

2.99 


3.87 







Concl usions 


During this contract Control Dynamics supported Marshall Space Flight 
Center in many capacities in promoting and achieving the goals involved with 
this project. The diversity of the tasks involved demonstrated the capabili- 
ties this company has available; the construction of an excellent structural 
model of the test structure (as seen when comparing the experimental and ana- 
lytical results), the review, analysis, and design of controllers to be used 
in controlling the test structure, the development of a simulation to use in 
runring test scenarios, and various miscellaneous tasks such as having the 
air bearing and tip bracket constructed. Control Dynamics hopes to be 
further involved with this innovative project in the future. 
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APPENDIX B. ATTACHMENT FROM MR 


# 13 PLANAR MODEL 


SIMPLIFIED PLANAR MODEL OF GROUND 


TEST VERIFICATION STRUCTURE 


Figure 1 depicts the planar version of the LSS-GTV c .tructure for which a 
model has been developed. The model is 12th order including 6 structural modes. 
One mode results from the base translation, one from the gimbal rotation, three 
from the flex characteristics of the Astromast Boom and gyro package, and one 
from the tether and tripod. The model Includes one Input, one disturbance and 
one output which are a control torque at the gimbal, a translational force at 
the base, and the rotation rate of the Astromast tip, respectively. 


The Gimbal flex program has been amended to output the system model In the 
following form 

n_ + + fi 2 n_ = dVj + cV£ * (1) 

0 = V3n . ( 2 ) 

The values n_ are the modal coordinates, d is the scalar disturbance input, c is 
the control input, and 0 is the tip rate sensor output. The values of omega, 
Vi, V2, and V3 along with the transition matrix from physical to modal coor- 
dinates (£ = Tn_, where £ are the physical coordinates) are attached. 


The model can then be cast into state variable form as 


ILL 


[ 0] 6 

*6 


HI 

+ 

06 

06 


d 

it 2 

u _ 


-H 2 

- 2 cfi 


H 2 


h 

V 2 _ 


c 


(3) 


D2J 

where nj = n and = n_. The terrn ^6 1s 3 6 x 6 zero matrix, 4 is a 6 x 6 
identify matrix, and O5 is a 6th order vector of zeros. By row and column 
interchanges, the system can be cast into block diagonal form so that the 
discretization process outlined in ED12-81-2 may be applied. Ongoing work 
includes production of coding to automate this process. 


0 = 


0 6 t V3 T 


B-l 
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BOOM 

(FIRST 3 MODES) 
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ATTACHMENT FROM MR # 14 ELIMINATION OF UNCONTROLLABLE RIGID BODY MODE 


ELIMINATION OF THE UNCONTROLLABLE RIGID BODY MODE 
FROM THE PLANAR SYSTEM MODEL 


The Gimbal flex program outputs the system model in the form 

+ 2cnn. + n 2 ji = V;d + 


( 1 ) 


and 



( 2 ) 


(See Progress Report 13.) The matrix, 0 , and vectors V^, ^ 2 , and V 3 , as given 
in the previous report, appear in Table 1 for convenience. The fifth and sixth 
diagonal elements of n reveal two very low frequency modes, relative to the 
ethers present; one at .0136 rad/sec and another at .0141 rad/sec. These are 
actually rigid body modes, however, for numerical reasons in computing, they do 
not appear as exactly zero frequency. 


Examination of Figure 1 reveal's that there are two rigid body modes modeled. 
One is due to the pure translation of the whole experiment assembly and is 
uncontrollable because the only control input is a torque at the gimbal. The 
other rigid body mode is due to the pointing of the gimbal and is the one we 
wish to control. The problem arises from the Vector V 2 which relates the 
control input to the modal coordinates. Since its fifth and sixth elements are 
both nonzero, neither rigid body mode can be thrown away as is. The two rigid 
body modes must be transformed so that one is unaffected by the control input 
and one is left controllable. The uncontrollable mode can then be eliminated 
from the model for design purposes. 


Consider a subspace of the modal coordinates containing only the rigid body 
modes (fijj = 0) and let it be denoted by the subscript "R". The system in 

this subspace is 

rjp = V^pd + V 2 R c (3) 

and 

T , . 

e R = V3RDR* 


Now the desired transformation is 

m = Rur- 


c-i 


(5) 


Premultiplying (3) by Q and then substituting according to (5) gives 


JQR = 9XlR d + 0y2R c • (6) 

For the case where there is exactly one controllable and one uncontrollable 
mode and a scalar Input, the requirement for Q is simply that QV£R have one 
zero element. Since there is some arbitrariness within this bound on the ele- 
ments of Q, let us require that Q be orthonormal so that we have the property 
Q' 1 = Q t and the Inverse is guaranteed to exist. 

To this end, let 


and 


911 

912 

921 

922 


V 2R 


b l 

b 2 


It is now obvious that 


1 b 2 

(bx 2 + b 2 2 )^ b! 

satisfies the requirements made. 


-bi 

b 2 


(7) 


( 8 ) 


(9) 


The output equation (4) may now be written in terms of the new modal coor- 
dinates, _np, and the transformation, Q, as 

0 = vJrQ-^. (10) 

Because of the special properties assumed for Q, this may be rewritten as 

0 = (QV3R) T 3iR (ID 

Now equations (6), (9), and (11) together describe the new rigid body subspace 
for the modal coordinate system so that one rigid body mode is not affected by 
the control input and one is controllable. 


Such a transformation was applied to the two rigid body modes of the original 
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planar system model and the results are shown In Table 2. Note that the order 
of the new coordinates was changed so that the rigid body mode appears as tyi 
which has been set to exactly zero. Some error is incurred in setting this 
value to zero; however, work is underway to alleviate the problem of nonzero 
"rigid body" modes being produced by the Gimbalflex program. 

Ongoing work Includes generalization of the process described here to any number 
of rigid body modes and control Inputs and automation of the process. 
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EARTH FIXED (LABORATORY FRAME) 


FIGURE 1. Planar Version of the LSS-GTV Structure. 
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TABLE 1. ORIGINAL 
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PLANAR SYSTEM MODEL 
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0.401554E-01 
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TABLE 2. PLANAR SYSTEM MODEL WITH UNCONTROLLABLE RIGID BODY MODE ELIMINATED 
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APPENDIX D. 

ATTACHMENT FROM MR #16, HP 9845C SPEED TEST RESULTS AND HARDWARE/SOFTWARE 


STUDY OF COMPUTATIONAL SPEED PROBLEM 



SPEED TEST ON HP 9845C 


The speed tests are summarized by Table 1 with the following definitions, 
n: system order 

m: number of outputs from digital controller 

r: number of inputs to digital controller 

The digital controller is implemented using matrix arithmetic and the 
following recurrence relation. 

x( i+1 ) = Ax ( i ) + Bu( i) 

y(i) = Cx( i ) 

i: iteration number 

A: nxn system matrix 

B: nxr input matrix 

C: mxn output matrix 

x: nxl state vector 

u: rxl input vector 

y: mxl output vector 
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TEST CASE* 


ms/ iteration 


Digital Controller (n=4, m=3, r=6) 14.5 

Digital Controller (n=9, M=3, r=6) 

vectors filled with zeroes- 46.0 

with write of values to floppy 47.0 

Digital Controller ( n=6, m=3, r=6) 

vectors filled with zeroes 27.7 

no zeroes in vectors 34.0 


♦all cases tested fc. 1000 iterations 


Table 1 . Summary of digital controller implementation 
speed test on HP 9845C. 
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As is obvious from Table 1, only the lowest order digital controller imple- 
mentation could be carried out within the required 20ms sample period. For this 
reason the HP 9845C desktop computer was deemed unfit for this particular task. 

Hardware/Software Study of Computational Speed Problem 

A study was made to determine the availability and approximate cost of com- 
puters and/or hardware accessories necessary to meet the 20ms sample period 
speed requirements. Additional requirements were that the control algorithym 
could be programmed in a high level language (FORTRAN OR BASIC) and that the 
machine have sufficient storage to store the data from a complete experiment. 

The results of the study are sumarized in Table 2 and informative comments 
with respect to the individual machines follow. 
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System w/30 I’.ASIC 

MP9020S | Mt> 5.25“ floppy lnter-n.il Internal GP-IO Graphics HP-UX 

(16-bit (single user) fUPTHAM 77 24 *As. J50.643 

S49.04S 10 Mb internal Monochrome Thermal parallel) w/Graphics/ 



Summary of liimln.irr/sof'.irarc requirements i „ meet c ompu I .• I mnal speed requirements for 
r. 1 V e * pe r l men I iliqilnl controller. (•Estimated cost) 









































INFORMATIVE COMMENTS 


1) The HP9020S is the single user, one package, desktop version of Hewlett 
Packard's new 32 bit processor machine. It includes hardware math and uses 
an 18MHz clock. Everything needed, hardware S software, with the exception 
of the GP-IO, is included with the ”S" package as shown in the table. 

2) The HP9020B is the same machine (package and processor) as the HP9020S with 
no standard peripherals or software. Shown in the table are the minimal 
add-ons required so that the system may meet the task requirements. 

3) The HPA700 (from the 1000 series) is a mini-computer geared toward multi- 
user environments and expandability. It is the lowest priced mini available 
from HP which takes advantage of hardware mathematics and should meet the 
speed requirements easily. 

4) The HP9845C system would be built around the existing 9845C in the test lab. 
The necessary purchase is essentially the array processor unit built by 
ANALOGIC. The array processor would speed up array manipulations only and 
would sup • no increase in speed for control algorithms which did not take 
advantage or' matrix arithmetic. There is much concern and uncertainty over 
the fact that in this system, the 9845 must communicate with both the COSMEC 
I and the array processor. This is seen as a possible deterrent to the 
increased speed gained by use of the array processor. 

5) Seattle Computer Company's Gazelle is an 8086 based micro system which takes 
advantage of the 8087 coprocessor. No benchmarks regarding speed are 
available from the company, however, considering the timing specs given by 
Intei for the 8086-8087 combination, the machine should be well within the 
speed requirements for the task. A dealer handling the Gazelle is located 
in Atlanta, and an inspection of the machine and possibly a test run are 
suggested. 


D-5 


| N8 5 13844 


SECTION 1 


SECTION 2 



APPENDIX E. 

METHOD OF DETERMINING MODE SHAPES AND NATURAL FREQUENCIES OF THE 
NASA UNMODIFIED TEST STRUCTURE 

CONTINUOUS BEAM CLOSED FROM SOLUTION TO THE NASA-LSS ASTROMAST 


TORSIONAL VIBRATION 



Appendix E 
Section 1 

Method of Determining Mode Shapes and Natural Frequencies of the NASA 
Unmodified Test Structure 

Abstract 

This paper illustrates the methods used to determine the lower natural frequen- 
cies and their corresponding mode shapes of the NASA-LSS Astromast (Unmodified 
Test Structure), and the mass integrals associated with the mode shapes. The 
test structure is modeled as a cantilever beam with 91 lumped masses and without 
the tip mass on the free end of the beam. This uncouples the torsion and 
bending modes and allows for them to be determined separately. We limited the 
frequency range to have an upper bound of 100. rad/sec (15.92 Hz.). In this 
range from 0.-100. rad/sec, we found three bending frequencies and one torsion 
frequency. 

Introduction • 

One of the simpler approaches available for determining eigenvalues and eigen- 
vectors is used in the development of this structural model. But use of a 
simplified method requires simplifying the assumptions being made. 
The first assumption is to treat the beam as cantilevered. This is a fairly 
safe assumption when considering the weight of the mast as compared to the 
weight of the LSS/GTY experiment equipment. The lumped mass approach is a 
second assumption which has been shown many times in the past to be reliable, 
especially when the distance between masses is relatively small as in this case. 
The ASTROMAST is modeled with no tip mass as this is how it is first tested in 
shaker tests. The tip mass is included in another program which incorporates 
the frequencies of the beam itself and uses modal synthesis to determine fre- 
quencies of the entire test structure. Without the tip mass/instrument package 
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connected to the free end of the beam, torsional and bending modes are 
uncoupled. This is due to there being no offset on the end which would couple 
the two modes. 


Problem Description and Method of Solution 

The problem in this case is to determine the natural frequencies, mode shapes, 
and mass integrals of the LSS Astromast. The beam is divided into 91 equal sec- 
tions, as in the actual mast itself, with mass equally distributed across the 
beam at each node and having essentially massless connectors between the nodes. 
Transfer matrices can now be generated to relate conditions from one point to 
another anywhere on the beam. The matrices for the bending and torsion cases 
are shown in Attachment 1. The procedure for determining eigenvalues in transla- 
tion is shown below. The torsional method is similar but much easier due to the 
smaller size of the transfer matrix. 

Starting with the tip (free end), the point and field matrices can be multiplied 
back to the base to yield the following equation: 


r h 

y 

1 



y 

e 


U1 


e 

M 




M 

V 




[yj 


T B 


Utilizing the fact that the moment and shear are zero at the tip and deflection 
and slope are zero at the base, a subsequent matrix can be formed relating known 
tip conditions to unknown base conditions: 



Therefore, to have a nontrivial solution of the above equation; the determinant 
of matrix U22 must equal zero. Since U22 is a function of the system's natural 
frequency, when the determinant is zero an eigenvalue has been located. 

Once an eigenvalue is defined, the mode shapes can be calculated. Going back to 
the above equation, set M-1.0 and then solve for V. This procedure generates a 
consistent set of base conditions for the structure. The transfer matrices can 
now be used again once the base conditions are established: 


y 


r— - 

'o' 

0 



0 

M 

= 

Tl 

1. 

V 



V 

V. - 



v» J 
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This equation solves for conditions at lump number 1 relative to the base. Then 
conditions at 2 can be solved based on conditions at 1, and this process is con- 
tinued until the tip conditions are reached. The intermediate values are writ- 
ten to d! o t files where they can later be retrieved and plotted as eigenvectors. 
In the same routine the mass integrals can also be calculated using the eigen- 
vector data and the following formulas for translation: 

To = Sum (mi*yi) 

Ti = Sum (mi*li*yi*) 

M = Sum (mi*(yi**2) ) . 

Basically, the same process is used for the torsional modes as that used for the 
bending modes. The tip and base conditions are now related by the following 
formula: 
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Only one term, which is a function of natural frequency, is now set to zero 
instead of a determinant of a matrix. The base torque is now set equal to 1.0 
and the values for twist and torque are determined from the individual transfer 
matrices just as in the bending case. The formulas used for the rotation mass 
integrals are as follows: 

To = 0. 

= Sum (Ii*<}>i) 

M = Sum (Ii*Ui**2)). 


Resul ts 

Using the methods described previously and a range of 0.-100. rad/sec, three 
bending modes and one torsional mode were located. The three bending modes 
occurred at 

3.882 rad/sec (00.62 Hz), 

24.609 rad/sec ( 3.92 Hz), 

70.238 rad/sec (11.18 Hz); 
and the torsional mode occurred at 
43.208 rad/sec ( 6.88 Hz). 

The plots for these mode shapes can be seen in Attachment 2. (Dr. George Doane's 
closed system analysis (Section 2) backed up the torsional frequency and an 
independent program written by Dr. Eugene Worley backed up the bending mode 
values. ) 
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The mass integrals calculated are as follows for each frequency: 


3.882 24.609 70.238 43.208 

r 0 1.725E-04 1 . 597 E -05 2.767E-06 0.000 

r : 6.481E-02 2. 128E-03 8.731E-05 5.265E-01 

M 3.779E-06 9.151E-08 1.069E-08 1.338E-06 

The bending mass integrals agreed very closely with Dr. Worley's calculated 
values. 

Conclusions 

Therefore, the methods and assumptions used in this paper are deemed feasible 
and very reliable and are recommended to be used in similar projects in the 
future. 
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Attachment 1. UNCOUPLED TRANSFER MATRICES 


Transl ation: 



Point 




Field 

~1 

0 0 

0 


1 

5.625 

Ay0 

0 

1 0 

0 


0 

1 

A0O 

0 2 

0 1 

0 


0 

0 

1 

mu 

0 0 

1 


0 

0 

0 


Rotation: 





p- 

Point 




Field 

[-Ju, 2 

9 



[j 

ATOR] j~<t> 





\y 

u. 



5.625 

1 




KTR 


1.80464 

5.07554 



-5.07554 • lol 
3.97252 • 10_ 


A = KTR- 
A( 1 ,1 ) = Ayy 

A(1 ,2) = A(2,l ) = Ay0 = AGy 
A(2,2) = A30 


ATOR = 5.58 • 10- 5 

Beam length = 512 inches 
Section length - 5.625 inches 
D = 4.457 inches 
Mass of beam = 0.012940 snails 

Mass of one lumped mass = (0 .012940 ) / ( 9 1 ^ = MI snails 2 
Inertia of one lumped mass = JI = MI • D snail s-inches 
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ATTACHMENT 2 .- MODE SHAPES 








v>\ LCJL'C'H v-Dcru <>IC PL MOHClUr 

o.oeoor+co o.coooo+oo o.oooon+oo o.ooooo+oo o. 100011+01 

O. 10000+01 0.36250+01 0.58070-05 0 . 25? 1 D -OS 0.90500 + 00 
0.20000 + 01 0.1 1250 + 02 0. 26030-0'! 0.51110-05 0. 96990 + 00 
0.30000+01 0.16990+02 0.60600-01 0.76570-05 0.95170+00 
0.10000 + 01 0.22500 + 02 C. 10910-03 0.10130-01 0. 92990 + 00 
0.50000+01 0.29120+02 0.17130-03 0.12560-01 0.92130+00 
0.60000+01 0 . 32750+02 0.21710-03 0.11960-01 0.90900+00 
0. 70000 +01 0. 39530 + 02 0,33670-03 0.17320-01 0. 87130 + 00 
0.80000+01 0.15000+02 0.13920-03 0.19620-01 0.87970+00 
0. 90000+01 0.50630+02 0.55150-03 0.21910-01 0.86170+00 
0.10000+02 0.56250+02 0.68260-03 0.21150-01 0.81970+00 
0.11000+02 0.61990+02 0.92310-03 0.26350-01 0.83170+00 
0.12000+02 0.67500+02 0.97590-02 0.28510-01 0.81970+00 
0 , 13000+02 0. 73130+02 0.11110-02 0.30630-01 0.80170+00 
0.11000+02 0.78750+02 0.13170-02 0.32710-01 0.78970+00 
0.15000+02 0.91370+02 0.15060-02 0.31750-01 0.77180+00 
0.16000+02 0.90000+02 0.17050-02 0.36760-01 0.75990+00 
0.17000+02 0.95630+02 0.19160-02 0.33720-01 0.71190+00 
0.18000+02 0.10130+03 0.21330-02 0.10650-01 0.73010+00 
0.19000 + 02 0.10690+03 0.23 700-02 0. 12530-01 0.71320 + 00 
0.20000+02 0.11250+03 0.26130-02 0. 11380-01 0.70010+00 
0.21000+02 0.11910+03 0.28670-02 0.16190-01 0.68550+00 
0. 22000+02 0.12380+03 0.31300-02 0.17960-01 0.67080+00 
0.23000+02 0.12910+03 0.31030-02 0.19690-01 0.63600+00 
0.21000+02 0.13500+03 0.36860-02 0.51390-01 0.61130+00 
0.23000+02 0.11060+03 0.39780-02 0.53010-01 0.62670+00 
0.26000+02 0.11630+03 0.12800-02 0.51660-01 0.61210+00 
0.27000+02 0.15190+03 0.15900-02 0.56210-01 0.59760+00 
0.28000+02 0.15750+03 0.19100-02 0.5778D-01 0.58310+00 
0.29000+02 0.16310+03 0.52370-02 0.59290-01 0.56860+00 
0.30000+02 0.16880+03 0.55710-02 0. 60750-01 0.55130+00 
0.31000+02 0. 17110+03 0.59180-02 0. 62190-01 0.53990+00 
0.32000+02 0.18000+03 0.62700-02 0.63570-01 0.52570+00 
0.33000+02 0.18360+03 0.66300-02 0.61920-01 0.51160+00 
0.31000+02 0.19130+03 0.69980-02 0.66210-01 0.19750+00 
0.35000+02 0.19690+03 0. 73720-02 0.67320-01 0.18350+00 
0.26000+02 0.20250+03 0.77510-02 0. 68770-01 0.16960+00 
0.37000+02 0.20810+03 0.81130-02 0.69990-01 0.15390+00 
0.38000+02 0.2138D+03 0.85390-02 0.71150-01 0.11210+00 
0.39000+02 0.21910+03 0.89110-02 0.72290-01 0.12830+00 
0 . 10000 + 02 0.22500 + 03 0. 93190-02 0.73390-01 0.11500 + 00 
0.11000+02 0.23060+03 0.97610-02 0.71130-01 0.10160+00 
0. 12000 + 02 0.23630 + 03 0.10180-01 0. 75180-01 0.38820 + 00 
0.13000+02 0.21190+03 0.10610-01 0. 76180-01 0.37310+00 
0.11000+02 0.21 75D+0Z 0.11010-01 0.77110-01 0.36210+00 
0.15000 + 02 0.25310 + 03 0.11130-01 0. 78370-01 0.31920 + 00 
0. 16000 + 02 0. 25880 + 03 C.11V2D-01 0.79270-01 0. 33650 + 00 
0.17000+02 0. 2611 D +03 0.12370-01 0.80130-01 0.32390+00 
0.18000+02 0.27000+03 0.12820-01 0.8096D-01 0 . 31110+00 
0.19000 + 02 0.27360+03 0.13280-01 0. 81750-01 0. 29910 + 00 
0.50000+02 0. 28120+03 0. 13710-01 0.82520-01 0.28700+00 
0.51000 + 02 0. 28690+03 0.1 1200-01 0. 83250-01 0.27500 + 00 
0.52000+02 0.29250+03 0.11670-01 0.83960-01 0.26320+00 
0.53000+02 0. 29910+03 0.15110-01 0.91A3D-01 0.25160+00 
0.51000+02 0.30380+03 0.15620-01 0.85270-01 0.21010+00 
0.55000+02 0.30910+03 0.16100-01 0.95980-01 0.22990+00 
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.OOPF02 
'000F02 
300D+C2 
9000 FC2 

0000+02 

• 1 0 0 0 F 0 2 
.2000+02 
>3000+02 
64000+02 
65000+02 
6 £ 0 0 0 F 0 2 
6 7 0 0 0 F 0 2 
.6800DF02 
•6V0OPF02 
•7000DF02 
•7100 UFO 2 
• 7200 OF 02 
i • 73000+02 
1.7400DF02 
0 7 3 0 0 0 F 0 2 
0.76000+02 
0.7700PF02 
0.78000+02 
0. 79000+02 
0.80000+02 
0 . 8 1 000+02 
0.8200D+02 
0 . 8 3 0 0 D F 0 2 
0 .84000+02 
0.83000+02 
0.8600D+02 
0 .8700PF02 
0 .89000+02 
0 . 8 V 0 0 LH 0 2 
0 . 90000+02 
0.V100DFC2 


0.3150DF03 
0 . 3 2 0 £ £' F 0 3 
0.32630+03 
0. 33190 F03 
0.33750+03 
0.3131 Li F03 
Q.34880+03 
0 . 3 3 1 1 Li F 0 3 
0 . 3 £ 0 0 1* F 0 3 
0.36560F03 
0 . 3712PF03 
0 .3769DF03 
0 . 3B25PF03 
0.3881 [i F03 
0.3938 Li F03 
0.39V4DF03 
0.4050PF03 
0. 1106DF03 
0 . 4163LH03 
0. 421VDF03 
0 . 42750F03 
0. 1331DF03 
0. 4388 0+03 
0. 41440+03 
0.45000+03 
0. 45560+03 
0.46130+03 
0 . 46690+03 
0.47250+03 
0.47810+03 
0 .48380+03 
0. 48940+03 
0.49500+03 
0 . 50060+03 
0.50630+03 
0.51200+03 


o . ldsen-oi 

0 . 1 707P-01 
0.17t)dU-01 
0.18050-01 
0.18550-01 
0.19050-01 
0. 1955D-01 
0.20050-01 
0.20560-01 
0. 21060-01 
C. 21570-01 
0 22080-01 
0.22600-01 
0.23110-01 
0.23630-01 
0.21110-01 
0 . 24660-01 
0 . 23180-01 
0.25700-01 
0.26220-01 
0.26740-01 
0.2/260-01 
0. 27790-01 
0.28310-01 
0.28840-01 
0. 29360-01 
0.29eBD-01 
0.30110-01 
0. 30940-01 
0.31160-01 
0.31990-01 
0.32510-01 
0.32040-01 
0.33560-01 
0.34090-01 
0.31620-01 


0.86470-04 
0 .87020-01 
0 . 87550-04 
0.89050-01 
0.88520-04 
0.83960-01 
0. 89380-04 
0 .99780-01 
0. 90150-04 
0.90190-01 
0.90820-04 
0.91 120-01 
0.91390-04 
0.91650-01 
0.91890-04 
0.92100-01 
0.92300-04 
0.92180-01 
0.92640-04 
0 . 9279D-0 1 
0.92920-04 
0.93030-01 
0.93130-04 
0.93220-01 
0 . 93290-04 
0.93350-01 
0 93410-04 
0.93450-04 
0.93480-04 
0.93510-01 
0 .93530-04 
0.93510-01 
0.93550-04 
0.93550-01 
0.93560-04 
0 .93360-01 


0.21760F00 
0 • 20690F00 
0 . 19620F00 
0 . 1858DF00 
0 .17560+00 
0 . 16560F00 
0 . 15560F00 
0. 1163DF00 
0. 13700 FOO 
0 . 12800F00 
0.1192DFC0 
0 . 1 106DFOO 
0 . 10240F00 
0. 91350-01 
0.86640-01 
0.7920D-01 
0.72060-01 
0.65210-01 
0 . 58670-01 
0.52410-01 
0.46520-01 
0. 10930-01 
0.35660-01 
0.30730-01 
0.26150-01 
0.21900-01 
0 .18020-01 
0. 11190-01 
0. 11330-01 
0.S511P-02 
0.61320-02 
0.11090-02 
0.24780-02 
0 .12450-02 
0 .41730-03 
0 . 15060-13 




THY Y (24.609 R/S) 



000 150*000 200.000 ' * 250.000 ' ~ SOOTOOO"" 350.000 

LENGTH (IN) 



FnE - GJ .- W'bltzCC 

MM • ULKi&TH \f 'bZ.VL. 

0 . 00000+00 0 . 0 0 0 0 1 1 + 0 0 


0 . 10000+01 
0 . 200011 + 01 
0.300011 + 01 
0 . >10000 + 01 
o . rjooon+oi 
0 .60000 + 01 
0 . 70000+01 

o.eoooo+oi 

0 . 90000+01 
0 . 10000+02 
0 . 11000+02 
0 . 12000+02 
0 .130011 + 02 
0 , 1 >1000 + 02 
0 . 15000+02 
0 . 16000+02 
0 . 17000+02 
0 . 18000+02 
0 . 19000+02 
0 . 20000+02 
0 . 21000+02 
0 . 22000+02 
0 . 23000+02 
0 . 2 > l001i + 0 2 
0 . 25000+02 
0 . 26000+02 
0 . 27000+02 
0 . 26000+02 
0 . 29000 +02 
0 . 30000+02 
0 . 31000+02 
0 , 32000+02 
0 . 33000+02 
0 . 3 + 000+02 
0 . 33000+02 
0 . 36000+02 
0 . 37000+02 
0 . 38000+02 
0 . 39000+02 
0. >10000 + 02 
0 . >11000 + 02 
0. >12000 + 02 
0 . >13000 + 02 
0 . ++ C0D + O2 
0 . >13000 + 02 
0. >16000 + 02 
0 . >17000 + 02 
0. >16000 + 02 
0 . >19000 + 02 
0 . 50000+02 
0 . 51000+02 
0 . 5200 D +02 
0 . 53000+02 
0 . 5 + 000+02 
0 . 33000+02 


0 . 36250+01 
0 . 11250+02 
0 . 168SD+02 
0. 2250H+02 
0 . 28120 + 02 
0 . 33750+02 
0 . 39380+02 
0 . + 5000+02 
0 . 50630+02 
0 . 56250+02 
0 . 61880+02 
0 . 67500+02 
0 , 73130+02 
0. 7875 0 + 02 
0 . 8 + 370+02 
0 . 90000+02 
0 . 93630+02 
0 .10130 + 03 
0 . 10690+03 
0 . 11250+03 
0 . 11810+03 
0 . 12380+03 
0 . 129 + 0+03 
0. 13500 + 03 
0 . 1 +060 + 03 
0 . 1 + 630+03 
0 .13190 + 03 
0 . 15750+03 
0 . 16310+03 
0 . 16880+03 
0 . 17 + + 0+03 
0 . 18000+03 
0 , 18560+03 
0 . 1913T+03 
0 . 19690+03 
0 . 20250+03 
C . 20810+03 
0 . 21380+03 
0 . 219 + 0+03 
0 . 22500+03 
0 . 23060+03 
0 . 23630+03 
0 . 2 + 190+03 
0 . 2 + 75D + C3 
0 . 25310+03 
0 . 25830+03 
0 . 26 ++ D+03 
0 . 27000+03 
0 . 27560+03 
0 . 28120+03 
0 . 28690+03 
0 . 29250+03 
0 . 29310+03 
0. 3 038 0 + 03 
0 . 309 + 0+03 


0 . 00000+00 
0 . 20 + 60-05 
0 . 1 8020 - 0 + 
0. +716 0 - 0 + 
0 . 88700 - 0 + 
0 . 1 + 190-03 
0 . 20600-03 
0 . 28020-03 
0 . 36380-03 
0 . + 5600-03 
0 . 5561D-03 
0 . 66330-03 
0 . 77690-03 
0 . 89620-03 
0 . 10200-02 
0 . 11 + 90-02 
0 . 12810-02 
0 . 1 + 150-02 
0 . 15520-02 
0 . 16900-02 
0 . 18290-02 
0 . 19630-02 
0 . 21060-02 
0 . 22 ++ 0-02 
0 . 23790-02 
0 . 25120-02 
0 . 26 + 10-02 
0 . 27670-02 
0 . 28890-02 
0 . 30060-02 
0 . 31180-02 
0 . 32250-02 
0 . 33250-02 
0 . 3 + 180-02 
0 . 350 + 0-02 
0 . 35830-02 
0 . 365 + 0-02 
0 , 37170-02 
0 . 37710-02 
0 . 38160-02 
0 . 38530-02 
0 . 38300-02 
0 . 3S97D-02 


5 > LOPE - 
0 . 00000+00 
0 . 25 + 30-05 
0 .+ 9520-05 
0 . 722 + 0-05 
0 . 93620-05 
0 . 11360-0 + 
0 . 13230 - 0 + 
0 . 1 + 960-0 + 
0 . 16560 - 0 + 
0 . 1 S 02 D - 0 + 
0 . 1935 D - 0 + 
0 . 20550-0 + 
0 . 21620 - 0 + 
0 . 22550 - 0 + 
0 . 23350 - 0 + 
0 . 2 + 030 - 0 + 
0 . 2 + 570 - 0 + 
0 .2 + 990-0 + 
C . 25290 - 0 + 
0 . 25 + 60 - 0 + 


MOME-lOT 
0 .10000 + 01 
0 . 9 + 320+00 
0 . 89630+00 
0 .8 + + 5 D + 00 
0 . 79270+00 
0 . 7 + 090+00 
0 . 68920+00 
0 .63760 + 00 
0 . 58620+00 
0 . 53 + 90+00 
0 .+ 8380+00 
0 . +3300 + 00 
0 . 38250+00 
0 . 332 + 0+00 
0 . 28280+00 
0 . 23360+00 
0 . 18500-100 
0 . 13700+00 
0 . 89690-01 
0 . -131 + D -01 


0 .25520-0+ -0.260 10-02 
0 .25 + 50-0 + -0. +7 + 50-01 
0 .25270-0+ -0.9 1350-01 
0 .2 + 9 70-0 + -0 .13 + 20 + 00 
0.2+ 570-04-0. 17600+00 
0.2 +06D-0+-0 .21660+00 
0 .23++ 0-0+-0 .25610+00 
0.2272D-0+-0. 29+20+00 
0. 2 1 91 D-0+ -0.33 100+00 
0 .20990-0 +-0. 366 +0+00 
0 . 1999D-0 + -0. + 00 + 11 + 00 
0. 1 89 10-0 +-0. +3280+00 
0. 1 77+ D-0+-0.+ 6370+00 
0 . 16+90-0 +-0 . +9290+00 
0 , 15160-0 +-0. 52050+00 
0 . 1377D-C+-0.5+6+D+00 
0 . 1 23 1D-0 + -0 .57060 + 00 
0 . 1 0790 -0 + -0 .59300 + 00 
0. 92 18H-05-0. 6 1360+00 
0 . 75910-05-0.632 + 0 + 00 
0. 59 18U-05-0 .6+930+00 
0 . +2030-05-0. 66 ++0+00 
0 .2+510-05-0.67760+00 
0.39050-02 0.66690-06-0.68830+00 
0. 3 903D-02-0.il +4 D-05-0. 69820+ 00 • 
0 .33900-02-0 .29770-05-0.70580 + 0 0 
0 .336 8 0-02-0 .+8270-05-0.711 + 0 + 00 
0. 33350-02-0. 6690 D-05-0. 71 5 10+00 
0. 37920-02-0 .85590-05-0 .71700 + 0 0 
0.3733 0-02-0 . 10 + 3D-0 + -0. 71710 + 00 
0.3 675 0 -02-0 . 1 230D-0 + -0 . 7 1 530 + 00 
0. 36010-0 2-0. 1 +1 70-0 +-0. 71 130+00 
0.35160-02-0. 1 602 D-0 + -0 .70650 + 00 
0.3+220-02-0 . 17350-0 + -0 .69950 + 00 
0.33170-02-0.19670-0+-0.6909D+00 
0 . 3202 0-02-0 . 21+6D-0 + -0.63C7D + 00 


Okie *- •- - 
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0.56000+02 0.31500+03 0 . 307BH-02-0 . 23220-04 -0 . 6687D+00 
0 . 57000+02 0 .32060+0 3 0 . 2V 4 4D-02-0 .24730-0 1-0 . 65560+00 
0.58000+02 0 . 32630+03 0 . 28001' -02-0 . 26 640- 04 -0 . 6 4 090 + 00 
0 .57000+02 0 .33190 + 03 0 . 26 470 -02-0 . 2829 D-0 1-0 . 62 13 D + 00 
0 . 60000+ 02 0.33750+03 0 . 24850-02-0 . 2990 Li-0'1 -0 . 60740+00 
0 .610011 +02 0.34310+03 0 . 23 1 4D-02-C . 31 17 D-0 1-0 . 53390+00 
0.62000+02 0.34880+03 0 . 21 35D-0 2-0 . 327811-0'! - 0 . 1567211 + 00 
0.43000 + 02 0.33440 + 03 0 . 1 9 43 D -0 2-0 . 3 4 4 4 D-0 4 -0 . 5 4 350+00 
0.64000+02 0.34000 + 03 0 . 17320-02-0 . 3584 D-04 -0 . 526SD+00 
0 . 43000+02 0.34340+03 0.13470-02-0.37190-04-0.50430+00 
0.44000+2 0.37130+03 0.13230-02-0.38470-04-0.48110+00 
0.47000 + 02 0.37470 + 03 0 . 1 1 2 1 D -02 -0 . 39 700-0 4 -0 . 45720 + 00 
0.43000+02 0.33250+03 0 . 8767D-03-0 . 4086D-04-0 . 43270+00 
0.47000+02 0 .33310+03 0.44420-03-0.41940-04-0.40770+00 
0.70000+02 0. 39380+03 0 . 42970-03-0 . 42990-04-0 . 38270+00 
0,71000+02 0.37740+03 0 . 18770-03-0 . 43940-04-0 . 35730+00 
0 . 72000 + 02 0.4050 0 + 03-0. 59550-04 -0.44 84 D-04 -0.331 70 + 00 
0.73000+02 0 .41040+03-0 .31170-03-0.45490-04-0.30430+00 
0.74000+02 0.41430+03-0.54840-03-0.44450-04-0.28100+00 
0.75000+02 0 .42190+03-0.32720-03-0 .47140-04-0.25570+00 
0.76000+02 0.42750+03-0.10940-02-0.47790-04-0.23130+00 
0.77000+02 0. 43310+02-0. 13620-02-0. 48360-04-0.20720+00 
0 .73000+ 02 0. 4 3SSD+ 03-0 . 16230-02-0. 48870-04-0 .16370+00 
0.79000+02 0.44440+03-0.19070-02-0.49320-04-0.16100+00 
0. 80000 + 02 0.4500 0 + 03-0. 21 S3 0-02-0. 49720-04-0 . 13930 + 00 
0.81000+02 0. 45540+03-0 ,24620-02-0.50050-04-0. 11 S60+00 
0.62000+02 0. 461 3D+03-0. 27420-02-0,50340-04-0. 99120-01 
0.83000+02 0. 46670 +03-0, 302 40-02-0. 505 7D-0 4 -C. 80960 -Cl 
0.84000+02 0. 47250+03-0. 33CSD-02-C. 50760-04-0. 64260-01 
0. 85000+02 0. 4 73 1D+03-0. 35720-02-0. 509 10-0 4-0. 49 170-01 
0.86000+02 0.48380+03-0.38780-02-0.51020-04-0.35820-01 
0.87000+02 0. 439 40+03-0 . 41640-02-0 .31100-04-0.24340-01 
0.68000+02 0.49500+03-0.44500-02-0.51150-04-0.14880-01 
C.87C0Dr02 0. 50060+03-0 . 473/0-02-0 . 51 1SD-0 4-0 .75310-02 
0.90000+02 0.50630+03-0.50250-02-0.51190-04-0.25730-02 
0.71000+02 0. 51 20D+03-0 .53130-02-0.51200-04 0.50530-14 
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LENGTH (IN) 



F£CQ.=10.£3b 

' LEtOGTVA V TbE-FL. “bUOPE. HOMCMT 
0 .0 GOOD +00 0. 00000+00 0.00000+00 0.00000+00 0.10000+01 
0.10000 + 01 0.56250+01-0.12230-05 0.250120-05 0.9161D + 00 
0. 2 000 D + 01 0.11250 + 02 0.1100U-01 0. 17810-05 0. 83220 + 00 

0 .20000 + 01 0.16380 + 02 0.35050-0-1 0.63-130-05 0.7-13-10 + 00 

0. -10000 + 01 0.22500 + 02 0.70920-0-1 0.66930-05 0.66-170 + 00 

0.50000+01 0.28120+02 0.11620-03 0. 10320-01 0.53120+00 

0.60000 + 01 0.33750 + 02 0.17000-03 0.11730-01 0.-19830 + 00 

0.70000 + 01 0.39390 + 02 0.23120-03 0.12920-01 0. -11600 + 00 

0. BO 00 0 + 01 0.-15000 + 02 0.29870-03 0.13900-01 0.33160 + 00 

0.90000 + 01 0.50630 + 02 0.37110-03 0. 11670-01 0.25110 + 00 

0. 10000+02 0.56250+02 0.11750-03 0.15230-01 0.17560+00 

0.11000+02 0.61380+02 0. 52660-03 0.15590-01 0.98670-01 

0.12000+02 0.67500+02 0.6073D-03 0.15750-01 0.23770-01 

0.13000 + 02 0.73130 + 02 0.6337D-03 0 . 1 572D-0 1-0 . 1873D-0 1 

0. 11000+02 0.78750+02 0.7696D-03 0 . 1 5500^01-0 . 1 1 850+00 
0.15000+02 0. 81370+02 0.91910-03 0 . 15100-01-0 . 13530+00 
0.16000+02 0 . 9000D+02 0. 92620-03 0.11510-01-0.21870+00 
0.17000+02 0.95630+02 0.1000D-02 0.13910-01-0.30310+00 
0. 18000+02 0.10130+03 0.10700-02 0 . 1 2930-01 -0 . 361 20+00 
0.19000 + 02 0.10690 + 03 0.11350-02 0 . 1 19 1 0 -0 1 -0 . 1 1 5SO + 00 
0 .20000+02 0.11250+03 0.11910-02 0.10770-01-0.16290+00 
0. 21000+02 0. 11810+03 0. 12170-02 0.95020-05-0.50530+00 
0.2200D+02 0.12380+03 0.12920-02 0 . 81 33D-05-0 . 51270+00 
0.23000+02 0.12910+03 0.13310-02 0.66710-05-0.57510+00 
0.21000+02 0.13500+03 0.13610-02 0 . 51370-05-0 . 60220+00 
0.25000+02 0.11060+03 0. 13810-02 0.35360-05-0.62100+00 
0.26000+02 0.11630+03 0.13970-02 0.18860-05-0.61020+00 
0. 27000+02 0.15190+03 0.11020-02 0 . 1 9990-06-0 . 65 1 00+00 
C. 28000+02 0.15750+03 0.13980-02-0.15070-05-0.65620+00 
0.29000+02 0.16310+03 0 . 1 3850-02-0 . 32200-05-0 . 65600+00 
0.30000+02 0.16880+03 0.13620-02-0.19250-05-0.65020+00 
0.31000+02 0. 17110+03 0 . 1 331 0-02-0 . 660SD-05-0 . 6391 0+00 
0 .32000+02 0.18000+03 0 . 1291 0-02-0 . 82560- 05-0 . 62270+00 
0.33000 + 02 0.18560 + 03 0 . 1 2 120-0 2-0 . 935 ID-05-0 . 60 1 30 + 00 
0. 31000+02 0. 19130+03 0 . 1 1 850-02-0 . 1 1 390-01-0 . 57190+00 
0.35000+02 0,19690+03 0 . 1 1 200-02-0 . 1 2S5D-0 1-0 . 51 390+00 
0. 3600 0+02 0. 20250+03 0 . 10170-02-0 . 1 1220-01-0 .50810+00 
0.37000+02 0.20810+03 0 . 967 1 0 -03-0 . 1 550D-0 1-0 . 16930+00 
0.38000+02 0.21380+03 0 . 88 120-03-0 . 1 6670-01-0 . 1 251 0+00 
0.39000+02 0.21910+03 0 . 7891 0-03-0 . 1 7720-0 1-0 . 37850+00 
0.10000 + 02 0. 22500 + 03 0 . 69 1 80-03-0 . 1 861 D-01 -0 . 32850 + 00 
0.11000+02 0.23060+03 0 . 53990 -0 3-0 . 1 9 130 -0 1 -0 . 27590 +00 
0.12000+02 0. 23630+03 0 . 1 8120-0 3-0 . 2008 D-01 -0 . 22070+00 
0.13000+02 0.21190103 0.37550-03-0.20530-01-0.16370+00 
0.11000+02 0.21750+03 0 .26150-03-0.20930-01-0.10530+00 
0.15000+02 0.25310+03 0 . 1521 0-03-0 . 2 1 130-0 1-0 . 15790-0 1 
0. 16000 + 02 0. 2583 0 + 03 0.39060-01-0.21170-01 0.11300-01 
0.17000+02 0.26110+03-0. 73750-01-0 .21030-01 0.71550-01 
0 . 1 800D+02 0.27000 + 03-0. 1S55D- 03-0. 20780-01 0.13150 + 00 
0.19000 + 02 0.27560 + 03-0. 29530-03-0. 20350-01 0.19370 + 00 
0.50000 + 02 0.2S12D + 03-0. 10250-03-0. 19770-01 0. 25180 + 00 
0.51000 + 02 0 . 286 9 0 + 03-0 . 50 60 0-03-0 . 19010-01 0.30330 + 00 
0. 52000 + 02 0.2925D + 03-0. 60530-03-0. 18160-01 0. 36270 + 00 
0.53000+02 0 .29810+03-0 . 69950-03-0 . 17130-01 0.11130+00 
0.51000+02 0.30380+03-0.78780-03-0.16000-01 0.16120+00 
0 . 35000+02 0.309 10 + 03-0. 86970-0. 4 - 0.11730-01 0*51010 + 00 

E “ 17 0f ilGlN Al p-. r - 


0.31500+03-0.91110-03-0.13310-01 
0.32060+03-0. 10 11D-02-0. 11010-01 
0. 3263 D +03-0. 10700-02-0. 10250-01 
0. 33 1 9 0 + 0 3-0. 11 20 0-02-0 .85720-05 
0.33750 + 03-0. 1161 Li - 02-0. 6819 0-05 
0.31310+03-0 . 1 1 V 2 D - 0 2 - 0 .‘30030-05 
0.31S8 0+03-0. 12130-02-0.31370-05 
0 .35110+03-0. 1221 U-02 
0.36000+03-0. 12250-02 
0.36560+03-0. 12150-02 
0. 371 30+ 03-0. 11950-02 
0 .376 V 0+03-0 . 1 1610-02 
0.38250+03-0. 1123D-02 
0.383 ID + 0 3-0. 10720-02 
0.39380+03-0.10110-02 
0. 3 V9 ID +03-0. 9 1020 -03 
0 . 10500+03-0 . 86000-03 
0 .11060+03-0.77070-03 
0.11630+03-0.67280-03 
0. 121 90+03-0. 56670-03 
0.12750+03-0.15300-03 
0 . 1331 D+0 3-0. 33220-03 
0.13S8D+03-0. 20190-03 
0 . 11110+03-0. 71520-01 
0. 67210-01 


0.56000+02 
0 . 57000+02 
0.58000+02 
0.59000+02 
0 .60000+02 
0 . 61000+02 
0 .62000+02 
0 .63000+02 
0 .61000+02 
0 . 65000+02 
0.66000+02 
0 .67000+02 
0 . 6800D+02 
0 .69000+02 
0,70000+02 
0 . 71000+02 
0 .72000+02 
0.73000+02 
0.71000+02 
0.75000+02 
0.76000+02 
0.77000+02 
0 .78000+02 
0.79000+02 
0 . 80000+02 
0. 81000+02 
0 .82000+02 
0.83000+02 
0.81000+02 
0.85000+02 

0 ,e6000+02 

0.87000+02 
0.88000+02 
0.89000+02 
0.90000+02 
0 .91000+02 


0.15000+03 
0.15560+03 
0 .16130+03 
0.16690+03 
0.17250+03 
0.17810+03 
0.18380+03 
0.18910+03 
0.19500+03 
0 *50060 + 03 
0.50630+03 
0.51200+03 


0. 21080-03 
0.35870-03 
0.51030-03 
0.66500-03 
0.82210-03 
0.98190-03 
0.11130-02 
0. 13060-02 
0.11700-02 
0.16350-02 
0.18000-02 


-0 . 12210-05 
0 .70030-06 
0.26130-03 
0 .15910-05 
0 . 652 ID-05 
0 . ei31D-05 
0 . 1031D-01 
0 . 1211D-01 
0 .13910-01 
0 .15610-01 
0 .17230-01 
0.16770-01 
0 . 20220-01 
0.21570-01 
0. 22810-01 
0 .23950-01 
0 .21970-01 
0 .25880-01 
0.26630-01 
0.27370-01 
0.27950-01 
0.28120-01 
0 .28300-01 
0. 29090-01 
0.29290-01 
0.29130-01 
0.29510-01 
0.29550-01 
0.29560-01 


OF POOrt QJALTi 


0.55320+00 
0. 59220+00 
0.62730+00 
0.65910+00 
0 . 68160+00 
0.70610+00 
0 . 72360+00 
0.73590+00 
0.71310+00 
0.71610+00 
0.71 100 +00 
0.73720+00 
0.72580+00 
0.71000+00 
0.6S99D+00 
0 . 66530 + 00 
0.63810+00 
0.60690+00 
0.57270+00 
0.53580+00 
0.19670+00 
0.15530+00 
0.11360+00 
0.37060+00 
0.32730+00 
0. 23130+00 
0 . 21210+00 
0.20130+00 
0.16260+00 
0.12610+00 
0. 93510-01 
0,61520-01 
0 . 10010-01 
0.20660-01 
0.71010-02 
0.11790-13 
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F1LLQ . - M3. ZO?) ^DlSCC 


fOM 

0 . 00000+00 
0 . 1 0000+01 
0.20000+01 
0.30000+01 
0 .40000+01 
0.5 000 Li + 01 
0 . 6000D+01 
0.7000 LI + 01 
0 . 8000D+01 

o.voooLi+oi 

0. 1000D+02 
0 . 1 1000+02 
0. 1200D+02 
0. 1300D+02 
0. 14000+02 
0 . 15000+02 
0 . 1 600D + 02 
0.17000+02 
0. 18000+02 
0. 1900 D +02 
0.2000 Li + 02 
0 . 2 1 0 0 D + 0 2 
0.2200D+02 
0.23000+02 
0.24000+02 
0 . 23000+02 
0.2600D+02 
0 . 2700D+02 
0.28000+02 
0.29000+02 
0.3000D+02 
0.31000+02 
0. 32000+02 
0.33000+02 
0. 34000+02 
0 • 33000+02 
0.36000+02 
0.37000+02 
0.38000+02 
0.39000+02 
0.40000+02 
0.41000+02 
0 .42000 + 02 
0 . 43000+02 
0 . 44000 + 02 
0. 4 3000 + 02 
0.46000+02 
0.47000+02 
0. 4S0CD + 02 
0. 4900D + 02 
0.50000+02 
0.51000+02 
0 . 52000+02 
0.53000+02 
0.54000+02 

a r* r ,\ i> • *i 


UElOGTlA 
0 . 00000+00 
0.36250+01 
0. 11250+02 
0 . 16880+02 
0.22500+02 
0.28120+02 
0.33750+02 
0.39380+02 
0 .45000+02 
0.50630+02 
0.56250+02 
0 . 61880+02 
0.67500+02 
0.73130+02 
0.78750+02 
0.84370+02 
0.90000+02 
0.93630+02 
0.10130+03 
0 .10690 + 03 
0.11250+03 
0.11810+03 
0.12380+03 
0.12940+03 
0.13500+03 
0.14060+03 
0.14630+03 
0.15190+03 
0.15750+03 
0.16310+03 
0. 168eD + 03 
0.17440+03 
0.18000+03 
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Definition of Terms 


NM 

LENGTH 

Y-DEFL 

SLOPE 

MOMENT 

PHI DEFLECTION 
TORQUE 


number of point mass along ASTROMAST 

location along ASTROMAST of the corresponding NM (inches) 

deflection of ASTROMAST at corresponding NM (inches) 

slope at corresponding NM (in/in) 

moment at corresponding NM (in-lb) 

angular rotation at corresponding NM (radians) 

torque at corresponding NM (in-lb) 
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Appendix E 
Section 2 

A Continuous Beam 
Vibration 


Closed Form Solution 


to the NASA-LSS Astroinast Torsional 


As one of the solutions to find the torsional vibration frequencies of the 
NASA-LSS beams the following continuous beam solution is offered. 

It is necessary to develop the equations of the continuous beam in order to 
understand the calculation of the numerical values required. Consider a beam 



i, 


TVr~riW 

l dt i 


-+ 


oT 


dx 


After Thompson, let x be measured along the length of the rod. The angle of 
twist in any length dx of the rod due to a torque T is 


dB= J- 
C 


dx 


whorp . i , 

d© = anale, 4ujist across ^rrnm 

X - 4orcy3e- applied +o 

C - CjonsVant of -prcrpor4iona\»^j 

Rearranging yields 

o>B_ __ JT_ 
at " C- 


S'' B _ | l 

c- ay 
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Thus the change in torque across the segment is given by 


^ i 
^ x 


6t 



dx 


This torque increment is balanced by the inertial reaction term associated with 
the segment. 


f 


dx 


J 


^d,c 


where J is the polar mass moment of inertia about the longitudinal axis of the 
rod, per unit length of the segment. 

Thus both C and J are per unit length values and the dx's may be divided out. 
Thus 

- C. ^ 6 4 -ut Sfe 

■ j ^>V X - "ST 1 ' 


To calculate C let 


o9 = -J- d x. 


T= C 


69 
d V 


now for a segment dx = 5.625", 


T - 

i 


C 

3. \oZ% 


d&. 


From Dr. Glaese's model, acro ss a segment the relationship 

t - i.i'iRB io M e- 


obtains: 



*5 .IrL^ 


- 1 Ti \ r i6 « IO- 


C- 100^6.58 


To determine J, recall that tne total moment of inertia associated with the 
segment is given by 


4 HSb 


mass - 


5.D 




-4 


- l.HZZ^IO sna'ils 


J segment = 1.422 x 10“** x (4 .46 ) 2 = 2.829 x 10~ 3 snail-inches 2 
therefore the J per unit length is 


J -- = 5.03 MO-^ 


Jz- 


* 10^ _ p 

5.0-b*IO' 4 “ C 


ooi t io s 


J 2 ,- l.sibt) t 10" 


Thompson states that the solution to 

-JC 1 -, 


5~e 


to be 


eCc?i - (P,*in(ai N ) - ■&.xi[ !j x))( C-iir. lot T)a»iwt) 
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The boundary conditions f or a uniform rod in torsional oscillation with one end 
fixed and the other end free are 


(.h u>ne.n \-D ) Q-O 

(£) ioWlo X - H , •Vorcyje.-O-, or 
^ B 




; "-0 


Because tne time origin is arbitrary the general solution may be taken to be 


9 - (A«b\n(^r)+ "B CDS n 


Clearly (1) results in B=0; (2) in requiring 


to* (»*)£) -o. 


This is satisfied when 


LOn^- — v ^>TT Sit 
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. 1 (n+ j n-O.i.i, 


Hence 
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co n -_ 1 C n+ ' c% S c b(n+--^ s ) 


Construct the table: 


n 

u>n 

4n 

O 

M3 Hg. 

C- /is. 

i 

i 

I30.C6 

2o. 'H 

£ 

2.11. 13 

34. 55 
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ATTACHMENT FROM MR # 22 RESULTS OF MODAL SYNTHESIS MODEL 



ATTACHMENT A. 


Results of Modal Synthesis Model 
Definition of nomenclature: 


beam tip: 

beam base: 

elevation gimbal : 
azimuth gimbal: 
roll gimbal: 
base: 


translates in X and Y directions 
rotates about X, Y, and Z axes 
translates in X and Y directions 
rotates abou^ X, Y, and Z axes 
rotates about X-axis 
rotates about Y-axis 
rotates about Z-axis 
translates in X and Y directions 
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modal frequency 0.0000000 rad/sec 


BASE TRANSLATION 


0.7208357 

0.0000000 

o.ooooooc 

ROLL GIMBAL ROTATION 


0.0000000 



AZIMUTH GIMBAL ROTATION 


0.0000000 



ELEVATION GIMBAL ROTATION 


0.0000000 



BEAM BASE TRANSLATION 


0. 7208357 

0.0000000 

O.OOOOOOC 

BEAM BASE ROTATION 


0.0000000 

0.0000000 

O.OOOOOOC 

TIP INSTRUMENT PACKAGE TRANSLATION 

0.7208357 

0.0000000 

o.ooooooc 

TIP INSTRUMENT PACKAGE ROTATION 


0.0000000 

0 . 0000000 

0 . ooooooc 

MODAL FREQUENCY 0.0000000 

RAD/SEC 



BASE TRANSLATION 


0.0000000 

0.7208357 

o.ooooooc 

ROLL GIMBAL ROTATION 


0.0000000 



AZIMUTH GIMBAL ROTATION 


0.0000000 



ELEVA1IUN GIMBAL ROTATION 


0.0000000 



BEAM BASE TRANSLATION 


0.0000000 

0.7208357 

o.ooooooc 

BEAM BASE ROTATION 


0.0000000 

0.0000000 

o.ooooooc 

TIP INSTRUMENT FACKAGE TRANSLATION 

0.0000000 

0 . 7208357 

o.ooooooc 

TIP INSTRUMENT PACKAGE ROTATION 


0.0000000 

0.0000000 

0.000000' 

MODAL FREQUENCY 0.0000000 

PAD/SEC 



BASE TRANSLATION 


0.0101980 

-0.0012932 

o.ooooooc 

ROLL GIMBAL ROTATION 


0.1107835 



AZIMUTH GIMBAL ROTATION 


0.0000000 



ELEVATION GIMBAL ROTATION 


0.0000000 



BEAM BASE TRANSLATION 


0.0101980 

-0.0012932 

o.ooooooc 

BEAM BASE ROTATION 


0.0000000 

0.0000000 

0.118783! 

TIP INSTRUMENT PACKAGE TRANSLATION 

-0.2215391 

0.0230927 

o.ooooooc 

TIP INSTRUMENT PACKAGE ROTATION 


0.0000000 

A. 0000000 

0.118783! 

MODAL FREQUENCY 0.0000000 

RAD/SEC 

. 


BASE TRANSLATION 


-0 • 2426753 

-0.0000825 

o.ooooooc 

ROLL GIMBAL ROTATION 


0.0075818 



AZIMUTH GIMBAL ROTATION 


0.0063638 



ELEVATION GIMBAL ROTATION 


0.0000000 



BEAM BASE TRANSLATION 


-0.0791644 

-0.0000825 

0.000000< 

BEAM BASE ROTATION 


0.0000000 

0.0063638 

0.0075811 

TIP INSTRUMENT PACKAGE TRANSLATION 

3.2468774 

0.0017931 

-0.001574! 

TIP INSTRUMENT PACKAGE ROTATION 


0.0000000 

0.0063638 

0.007581! 

MODAL FREQUENCY 0.0000000 

RAD/SEC 



BASE TRANSLATION 


-0.0000439 

0.2133059 

0.000000' 

ROLL GIMBAL ROTATION 


0.0009661 



AZIMUTH GIMBAL ROTATION 


0.0000033 



ELEVATION GIMBAL ROTATION 


0.0064128 



BEAM BASE TRANSLATION 


0.0000413 

0.102556? 

0.000000' 

BEAM BASE ROTATION 


0.0064128 

0.0000033 

0.000966 

TIP INSTRUMENT PACKAGE TRANSLATION 

-0.0001019 

-3.2637730 

0.012510' 

TIP INSTRUMENT PACKAGE ROTATION 


0 . 0064 128 

0.0000033 

0.000966 

MUDAL FREQUENCY 5.913205 

RAD/SEC 



BASE TRANSLATION 


0.6340323 

0.0000961 

0 . 000000 ' 

ROLL GIMBAL ROTATION 


0.0017399 



AZIMUTH GIMBAL ROTATION 


-0.0497098 



ELEVATION GIMBAL ROTATION 


0 . 0000120 
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BEAM BASE TRANSLATION 


-0.6432107 

-0.0001104 

0 . 000000 

BEAM BASE ROTATION 

TIP INSTRUMENT PACKAGE IRAN SLAT 

TIP INSTRUMENT PACKAGE ROTATION 

ION 

0.0000120 
0.5985857 
-0 . 0000228 

-0.049709B 

-0.0000565 

0.0345617 

0.001739 

-0.008594' 

-0.016302 

MODAL FREQUENCY 7. 293426 

RAD/SEC 



BASE TRANSLATION 


-0.0003125 

0.4891170 

0 .000000 

ROLL GIMBAL ROTATION 


-0.0001548 



AZIMUTH GIMBAL ROTATION 


0.0000250 



ELEVATION GIMBAL ROTATION 


0.0602246 


«• 

BEAM FACE TRANSLATION 


0.0003297 

-0.5509624 

0.000000 

BEAM BASE ROTATION 


0.0602246 

0.0000250 

-0.000154 

TIP INSTRUMENT PACKAGE TRANSLATION 

-0.0003538 

0.5674684 

-0.088492 

TIP INSTRUMENT PACKAGE ROTATION 

• 

-0.0453606 

-0.0000090 

0.000052 

MODAL FREQUENCY 11.29450 

RAD/SEC 



BASE TRANSLATION 


-0.0302741 

-0.0009661 

0 . 000000 

ROLL GIMBAL ROTATION 


0.0184101 



AZIMUTH GIMBAL ROTATION 


0.0024533 



ELEVATION GIMBAL ROTATION 


-0.0000814 



BEAM BASE TRANSLATION 


0.0327604 

0.0004400 

0 . 000000 

BEAN BASE ROTATION 


-0.0000814 

0.0024533 

0.018410 

TIP INSTRUMENT PACKAGE TRANSLATION 

0.0370206 

-0.0081332 

0.000938 

TIP INSTRUMENT PACKAGE NUTATION 


. -0.0000700 

-0.0043433 

-0.715111 

7 MODAL FREQUENCY 20.44361 

RAD/SEC 



BASE TRANSLATION 


-0 • 40591 15 

-0.0021554 

0.000000 

ROLL GIMBAL ROTATION 


-0.0001958 



AZIMUTH GIMBAL ROTATION 


0.0339763 



ELEVATION GIMBAL ROTATION 


-0.0002847 



BEAM BASE TRANSLATION 


0.4670748 

0.0027617 

0.000000 

BEAM BASE ROTATION 


-0.0002847 

0.0339763 

-0.000195 

TIP INSTRUMENT PACKAGE TRANSLATION 

0.3724887 

0.0012984 

-0.023614 

TIP INSTRUMENT PACKAGE ROTATION 


-0.0003839 

0.0924256 

0.030584 

HUDAL FREQUENCY 23.29484 

RAD/SEC 



BASE TRANSLATION 


0.0019001 

-0.4684157 

0.000000 

ROLL GIMBAL ROTATION 


0.0001902 



AZIMUTH GIMBAL ROTATION 


-0.0001598 



ELEVATION GIMBAL ROTATION 


-0.0616555 



BEAM BASE TRANSLATION 


-0.00220 

0.5963746 

0.000000 

BEAM BASE ROTATION 


-0 . 04x6555 

-0.0001598 

0.000190 

TIP INSTRUMENT PACKAGE TRANSLATION 

-0.0018192 

0.3291788 

-0.179170 

TIP INSTRUMENT PACKAGE ROTATION 


-0.0919095 

-0.0005578 

-0.009847 

MODAL FREQUENCY 57.86488 

RAD/SEC 



BASE TRANSLATION 


-0.1169710 

-0.0530800 

0.000000 

ROLL GIMBAL ROTATION 


0.0000152 



AZIMUTH GIMBAL ROTATION 


0.0105723 



ELEVATION GIMBAL ROTATION 


-0.0074173 



BEAM BASE TRANSLATION 


0.1546738 

0.0750166 

0.000000 

BEAM BASE ROTATION 


-0.0074173 

0.0105723 

0.000015 

TIP INSTRUMENT PACKAGE TRANSLATION 

-0.1921383 

-0.0510709 

0 .102847 

TIP INSTRUMENT PACKAGE ROTATION 


0.0371458 

-0.122797 7 

-0.002712 

MODAL FREQUENCY 58.02358 

RAD/SEC 



BASE TRANSLATION 


0.0377894 

-0.1626017 

0 . 000000 


F-3 



ROLL ti TUBAL ROTATION 
AZIMUTH GIHBAL ROTATION 
ELEVATION OIHBAL ROTATION 
REAM RASE TRANSLATION 
REAM BASE ROTATION 

TIN INSTRUMENT PACKAGE TRANSLATION 
TIP INSTRUMENT PACKAGE ROTATION 


0.0000324 
-0.0034170 
-0.022/311 
-0 . 05000B1 
-0.0227311 
0.0637058 
0 . 1153053 


0.22VV638 

-0.0034170 

-0.161353V 

0.0400853 


0 . 000000 - 

0.000052 

0.213034 

-0.001361 


TRANSLATIONS (INCHES) 
ROTATIONS (RADIANS) 
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APPENDIX G. AIAA PAPER 

DEFINITION OF GROUND TEST FOR LARGE SPACE STRUCTURE (LSS) 


CONTROL VERIFICATION 



I. INTRODUCTION 


As the United States moves into the Shuttle era of space technology, there 
are i.jmerous proposals from the scientific, civilian and Defense Communities 
which envision the use of Large Space Structures (LSS ) . By definition, a LSS is 
very flexible, probably lightly damped, and exhibits multiple vibrational modes 
of very low frequency. 

Many of the missions alluded to above require high performance from the 
LSS. This high performance takes the form of such things as extremely accurate 
pointing of lenses and the attainment of vibration free observation image 
planes. 

To meet the rigorous control requirements of LSS, a new or extended body of 
control theory has been developed. To test these various schemes, MSFC (Marshall 
Space Flight Center) is establishing a LSS laboratory in which experimentation 
with large beams, LSS components, or even perhaps a given full-size LSS may be 
performed. 

This paper delineates the program currently underway to develop the labora- 
tory facility and concurrently to develop and install the first experiment in 
the laboratory. 


G-l 



II. LSS Ground Verification Experiment 


The Ground Test Verification (GTV ) experiment is described by the drawing 
of Figure 1. The first test article will be tne ASTROMAST beam as shown. The 
ASTROMAST is extremely lightweight (about 5 pounds) and approximately 45 feet in 
length and is constructed almost entirely of S-GLASS. It is of the type co be 
flown in SAFE I. 

The test article will be mounted on the faceplate of the Advanced Gimbal 
System (AGS) engineering model which, along with an additional torque actuator 
in azimuth, provides the control inputs for the system. The azimuth gimbal also 
provides a means of rotating the entire experiment manually to produce different 
test scenarios. The ASTROMAST will be gravity unloaded by a constant tension 
cable whose upper end is free to translate. 

The AGS will sit atop the base of the test structure which will be sup-, 
ported by air bearings or oil film bearings and will include actuators to pro-, 
vide translational disturbance inputs to the test fixture. These disturbances 
could represent astronaut pushoff or RCS (Reaction Control System) thruster 
fi ring. 

Six separately packaged inertial measurement assemblies comprise the 
control system sensors. Two of the packages, containing three axis translational 
accelerometers, are identical. One will be mounted on the mast tip, and the 
other on the base of the test fixture. Three other packages contain ATM (Apollo 
Telescope Mount) rate gyros and will be installed on the AGS faceplate. The 
sixth package, the Kearfott Attitude Reference System (KARS), will oe placed at 
the mast tip along with the accelerometer package. 



The signals from these instruments will be read by the C0SME«' I data 
gathering and control system, and processed according to the control strategy 
under scrutiny. The control actuacor signals will then be transmitted to the AGS 
as inputs to the dynamical system. 

The COSMEC I will be interfaced to a Hewlett Packard HP9845C desktop com- 
puter which will store data as it is collected during a test run, and then pro- 
vide post experiment dat* reduction and display offline. The controller inputs 
and outputs (measurements and commands) will be recorded at each sample period 
or at some multiple of sample periods. 
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III. Project Management 


The project management encompasses the technical and business management 
activities. These activities are required to plan, execute, control and report 
technical performance of the LSS ground verification experiment. Two more 
important and time consuming efforts, relative to the LSS ground experiment 
management, are developing and coordinating the schedules among the various par- 
ticipating HSFC elements and utilizing the resources in a cost-effective manner. 
To this end, the project management has been partitioned into three categories. 
They are planning and control, configuration management, and procurement 
management. 

The planning and control element provides for the planning, authorizing, 
and controlling of the LSS ground test effort. This activity also provides for 
timely visibility into the performance, cost, and schedule and is keyed to the 
appropriate f-SFC organizations and to the Work Breakdown Structure (WSS) and tr.e 
WBS dictionary. Integrated project schedules for the overall LSS are included 
as a function of this element. The control function integrates cost, schedule, 
and performance and relates progress and variance from the initial planning. 

The configuration management category provides for the definition and 
implementation of configuration management activities for both hardware and 
software required in the performance of the LSS experiment. In addition, it 
includes as-designed and as-built configuration (for all levels), change 
control, change tracking, and configuration accounting. 
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The procurement management element includes responsibility for providing 
the project performance survei l lance, cost control, and status reporting on all 
elements. It' also includes arranging for the acquisition of the test article 
and various components of the experiment. 

The integration of these elements into a smooth operation is the main 
objective of the project management. Once the proper integration is established 
in the various MSFC participating elements, constant follow-up with par- 
ticipating elements helps ferret out any particular problem areas. This proce- 
dure provides an excellent method for interaction and tor problem solving in 
the project management. 
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IV. Systems Engineering and Integration 


The scenario for the first experiment in the ISS/GTV experiment facility 
involves the ASTROMAST as the test article atop the AGS and employs the COSMEC I 
system as the digital controller. From a systems point of view, the entire 
apparatus; including the base, AGS, ASTROMAST, tripod, and measurement devices; 
is considered the plant, or process, and the COSMEC I and its associated hard- 
ware and software are the controller. 

Because the experiment employs a variety of types of measurement devices 
the COSMEC I system must deal with several different signal types as is 
expressed by the captions in Figure 1. To enhance the organization and speed of 
operation of the system, the measurement devices are monitored by hardware cards 
peripheral to the COSMEC I processor itself. Each device is sampled at a rate of 
50 samples/second, thus ensuring adequate bandwidth of the digital control 1 er./ 
The digital controller is discussed in greater detail in the section de> icattd ■. 
to the COSMEC I subsystem. 

Command outputs from the control algorithm Implemented in the COSMEC I are X 
converted to analog signals which drive the gimbal torquer amps. This conversion 
is handled by hardware cards analogous to tho r .e which monitor the measurement 
devi ces . 

In addition to the candidate control algorithm, the COSMEC I carries out an 
inertial strapdown algorithm for the measurement devices which takes into 
account the effects of the acceleration due to gravity and the angular displace- 
ment due to earth rotation. This algorithm conditions the measurements so that 
they are true with respect to the laboratory reference frame. 
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A typical experiment would include disturbing the system via one of the 
actuators at the base and measuring the results with the system operating open 
loop, i.e., with the digital controller not included, then disturbing the system 
in a similar manner but with a candidate control algorithm in place. A com- 
parison of the recorded results should give an indication of the success of the 
control algorithm. More importantly, a series of such tests could be used to 
determine the relative merit of different control schemes. 


V. Design, Development, Fabrication, and Checkout of Subsystems 

The subsystems which comprise the LSS/GTV experiment fixture as described 
are currently in various stages of development at NASA MSFC. Thorough verifica- 
tion of each of the subsystems in controlled test environments comprises a 
significant part of the preliminary system testing. The subsystems will not only 
be tested individually, but will be tested in an integrated laboratory environ- 
ment where each of the subsystems will interact with the others in a manner much 
as if an actual test article were being used. Such testing is designed to ensure 
proper operation cf the complete test fixture upon assembly. 

A. AGS Gimbal Drive 

The Advanced Gimbal System (AGS) is a precision, two axis gimbal system 

designed for high accuracy pointing applications. The AGS gimbals serve the 

elevation plane and a third gimbal has been added to the system in the azimuth. 

(See Figure 2.) The AGS essentially rides on the azimuth gimbal so that the AGS 

) 

and its performance remain unaltered by the presence of the third gimbal. The 
AGS receives torque commands from the COSMEC I data and control system in the 
form of analog inputs over the range of -10 to +10 volts. This saturation repre- 
sents the current limit of 27 amps which is built into the AGS servo amplifier 
as a protective measure. Because the AGS servo amplifier outputs a current 
which causes an applied torque proportional to the current, the control 
algorithms used in the COSMEC I must be designed to produce torque command 
signals. 

The AGS gimbal torquers, with the power supply and servo amplifiers used in 
the ground test verification experiment, can generate 37.5 ft. -lbs. of torque 
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over an angular range of approximately ± 30 degrees. The azimuth torquer is 
capable of generating 13.8 ft. -lbs. over an angular range of about ± 5 degrees. 
It can, however, be set manually to allow the ± 5 degrees of rotation at any 
position about the 360 degrees of azimuth freedom. This allows the test article 
to be rotated to any position desired without remounting. 

The power supply for the gimbal system (pitch, yaw, and roll) is a 100 amp 
d.c. supply with output voltage variable from zero to 35 volts. (The AGS 
requires 30 volts.) The line to the power supply is 203 volts, three phase, and 
in the worst case the input power to the supply should net exceed 2550 watts. 


B. C0SMEC I Data and Control System 

1. System Aspects 

The COSMEC I is an AIM 65 based micro system which is used to handle data' 
from the control system sensors, output commands to the control system 
actuators, transmit data for storage 'to the HP 9845C desktop computer, and 
implement the control and inertial strapdown algorithms. The COSMEC I is a 
powerful eight bit micro system which uses special hardware and software to 
allow the handling of a variety of devices (sensors, actuators, etc.) in real 
time. This, along with the use of high speed hardware arithmetic processors to 
reduce compu*v.ation time, makes the COSMEC I an excellent machine for use as a 
digital controller. 

2. Hardware 

The AIM 65, which is the heart of the COSMEC I, employs the Motorola MC6502 
mi coprocessor operated at a two megahertz clock rate. (See Figure 3.) The 
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system includes 32 kbytes of random access memory (RAM), an alphanumeric 
keyboard, a single line display, a cassette tape machine for mass storage, and a 
small printer. The entire system is housed in a very portable package much 
resembling a suitcase. 

The COSMEC I "reads" a variety of types of sensor output signals via inter- 
face cards which are an integral part of the COSMEC I system. (See Figure 4.) 
These cards allow the COSMEC I processor to interface in a similar manner (with 
regard to 'computation) with the ATM (Apollo Telescope Mount) rate gyros, the 
KARS (Kearfott Attitude Reference System), the accelerometer packages, and the 
AGS (Advanced Gimbal System), each of which has a different type output or input 
signal. The COSMEC I also features a real time clock which will prove useful in 
the recording of experimental data. 

In order to carry out the large number of calculations required for imple- 
mentation of the inertial strapdown algorithm and the. control algorithm, the 
COSMEC I employs four hardware arithmetic processors connected on the system 
bus. Each of these processors can execute a 32 bit floating point multiply in 42 
microseconds and they are operated so that they process in parallel, thus mini- 
mi zing computation time. The dynamic range of the processors is ± 9.2x10 so 
there is little possibility of exceeding the computational range of the machine. 
Also, this eliminates the need for scaling of measurements in order to avoid 
machine overflow. Using the arithmetic processor units and assembly language 
programming, the inertial strapdown algorithm can be executed in approximately 
ten milliseconds and the first proposed control algo. .thm in about six millisec- 
onds. This puts the total computation time at well under the allowed 20 
milliseconds required to meet the 50 hertz sample rate. 
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3. Software 

The software used in the COSMEC I system may be separated into fc-r basic 
groups: 1) utility software for handling the various hardware cards which inter- 
face to instruments, 2) software to implement the control algorithm, 3) software 
to implement the inertial strapdown algorithm, and 4) initialization and startup 
software to ready the instruments and equipment for a test. 

The hardware cards which interface the COSMEC I's processor to the measure- 
ment instruments and actuators are individual by their very nature, and some 
special software is required to handle each card. However, each card makes 
information available to the processor as digital words, which is the unifying 
feature of the system. 

The digital controller software for the first ground test experiment will 
implement a linear discrete multi variable controller having multiple inputs and' 
outputs. The controller will be in state variable form and will be programmed so 
that the system matrices are initial input data' to the program and can be stored 
on tape and easily changed. The first controller software will be designed to 
implement a controller of up to ninth order having nine inputs and three out- 
puts. 


Because the inertial measurement instruments measure with respect to iner- 
tial reference space, there is a natural bias in the measurements due to the 
acceleration of gravity and earth's rotation. That is, in the earth based 
experiment the accelerometers measure about one g acceleration downward and the 
rate gyros measure about 15 deg/hr rotation while at rest with respect to the 
laboratory reference frame. The inertial strapdown algorithm provides a means of 
removing this bias from the measurement instruments. 
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In order to give the measurement instruments initial conditions and begin 
measurement for a test, initialization software is provided for the instrument 
strapdown algorithm. To begin a test, the structure will be stabilized with 
respect to the laboratory reference frame and the initialization routine will be 
executed. The strapdown algorithm will then be started and the apparatus will 
be ready to carry out a test. 

C. Inertial Measurement Assemblies 

Three different types of inertiai measurement assemblies are planned for 
use on the Ground Test Verification structure in the first experiment. The 
Kearfott Attitude Reference System (KARS), the Apollo Telescope Mount (ATM) rate 
gyros, and two acccelerometer packages developed by NASA. Each of the instrument 
packages generates signals in a particular form different from the other instru- 
ments as was mentioned in the section dealing with the COSMEC I interface cards. 
These different signal types are pursued as each instrument package is discussed 
in detail in the following paragraphs. 

1. Kearfott Attitude Reference System 

The Kearfott Attitude Reference System (KARS) is an attitude measurement 
system designed for use in the U.S. Army remotely piloted vehicle. (See Figures 
5 and 6.) It provides measurement resolution of 13.9 X 10~3 deg/sec in the 
pitch and yaw axes (axes transverse to the ASTROMAST) and 25.0 X 10~3 deg/sec 
in the roll axis (axis along the length of the ASTROMAST). The dynamic range of 
the rate gyro outputs of the KARS is 40 deg/sec in pitch and yaw and 70 deg/sec 
in roll. Because of its light weight (8.9 pounds), the KARS will be used as the 
mast tip rotation sensor in the first ground test experiment. Although the KARS 
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includes accelerometers and outputs measurements of linear acceleration, the 
measurements are not used in the ground test experiment because of inappropriate 
scaling of the instruments. 

The output signals of the KARS are in the form of asynchronous digital 
pulses. One signal, the change in angular position in yaw for instance, requires 
two channels; one for pulses representing positive rotation and the other for 
pulses representing negative rotation. The-COSMEC I system accumulates the 
pulses over a 20 millisecond period to produce measurements of the angular rate 
and position of the ASTROMAST tip. 

The KARS outputs three digital health check signals. One signal represents 
a check of the motor voltage and the pickoff excitation voltage and the other 
two are over and under temperature signals. In the ground test experiment, these 
signals will be monitored at the system console. 

2. ATM Rate Gyros 

The Apollo Telescope Mount (ATM) rate gyro packages are designed to measure 
small angular rates very precisely. Each package measures angular rate in. one 
axis (See Figure 7.) with resolution finer than 0.5 X 10 _ 3 deg/sec and offers a 
dynamic range of ± 1.0 deg/sec. The ATM rate gyro packages will be mounted on 
the faceplate of the engineering AGS so that they will measure the rotation of 
the base of the test article. Because they will be rigidly mounted to the gimbal 
system, these sensor packages' wi 1 1 most nearly measure the actuator inputs to 
the system. 
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The output signals of the ATM rate gyro packages are ± 45 volt analogs and 
are handled by the analog to digital converter card of the COSMEC I system where 
they are converted to 12 bit binary words. 

The ATM rate gyro packages require a warmup period of about 20 minutes and 
then stabilize about 20 minutes after the warmup period ends, thus requ : ring on 
the order of 40 minutes to become ready for a test. Each package requires 1.5 
amps during warmup and then 1.25 amps after stabilization; both at 28 Vdc. 

3. Accelerometers 

Two identical accelerometer packages (See Figure 8.) will be used on the 
ground test experiment fixture in the first test. One package will be placed on 
the mast tip along with the KARS, and the other on the test fixture base as shown 
in Figure 1. The necessary electronics for each accelerometer package is 

included on board the instrument package itself as shown in Figure 9. 

* 

% 

• The accelerometers provide resolution finer than 0.0001 g and , a dynamic 
range of ± 3 g with a bandwidth of 25 to 30 Hz. They require about 20 minutes 
for warmup, during which time each package requires 1.2 amps at 28Vdc. After 
warmup the power requirement reduces to about .9 amp per package. 

The signals from the accelerometers are different from either those of the 
KARS or the ATM rate gyros. As in the case of the KARS, two channels are 
required for each of the degrees of freedom of the accelerometer package, i.e., 
six channels per accelerometer package. One channel of each pair carries a 2.4 
kHz square wave synchronization signal and the other channel carries the acce- 
leration information. Zero accleration is represented by a signal identical to 
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that of the synchronization channel, positive acceleration by an increase in 
frequency, and negative acceleration by a decrease in frequency as compared to 
the synchronization channel. As in the cases of the other instruments, these 
signals are monitored by a hardware card in the COSMEC I system. 


D. Beam Containment Structure 

The beam containment structure includes the tip support mechanism, the 
base, the disturbance actuators and signal sources, and power supply for the 
entire test fixture. This is essentially that equipment required of a laboratory 
to carry out dynamic testing of structures such as the ground test experiment. 
As depicted in Figure 1 the beam containment area can accommodate structures 
approaching 120 feet in height. Access is provided at various levels along the 
structure via catwalks. Also the control room for experimental operations is at 
a position about 50 feet above the base cf the structure., thus making the 
experiment quite accessible with regard to viewing the experimental operations 
and so forth. 

1. Tip Support Mechanism 

The tip support mechanism is a tripod mounted on air bearings as shown in 
Figure 1 with a constant tension cable extending downward to connect to the tip 
instrument package. The cable connection point is incorporated into the tip 
instrument package mount assembly. The tip support mechanism is designed to 
support the weight of the tip instrument package and possibly some of the weight 
of the ASTROMAST beam without affecting the structural properties of the system 
by adding stiffness to the beam tip motion. The tripod is allowed about 13 
inches of total translational freedom in each of its two translational axes. 
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This will prove to be quite adequate for the GTV experiment because the tip 
deflections are expected to be less then six inches and the distance between the 
mast tip and the tripod support floor is about 60 feet as compared to the rela- 
tively short tripod of three feet. 

2. Base 

The base for the GTV experiment will be supported by air bearings in a 
similar manner to the tip support tripod. It is designed to be restrained as 
little as possible in the two translational degrees of freedom but not allow any 
rotational freedom. Translation of the base provides the means of applying 
disturbances to the system. The disturbance actuators will be hydraulic and/or 
electromechanical . in nature depending upon the frequency band of the applied 
disturbance signal. The disturbance system is designed to allow independent 
input signals to the two translational axes whether they be broadband random, 
sinusoidal, or a combination of both. 


E. Facility Computer Interface and Utility Software 

During testing, the COSMEC I data and control system will send data for 
storage to the HP 9845C desktop computer via a 16 bit parallel interface. Once a 
test run is completed, the HP 9845C will be available for display of data and 
possibly some data reduction. Also available is a Hewlett Packard 1000 series 
minicomputer to which the data stored on the HP 9845C can be transferred. The HP 
1000 series mini system includes the necessary software for several sophisti- 
cated data reduction techniques. Using the combined power of the two machines, 
much will be possible in the way of data reduction and display for post analysis 
efforts. 
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VI. System Model and Expected Performance 

For purposes of system studies and controller design, an analytical model 
of the Ground Test Experiment is necessary. This section describes the develop- 
ment of such a model and the following section compares it to tie results of 
modal tests of the ASTROMAST beam carried out at NASA MSFC. 

A. Introduction 

Modeling of the Ground Test Verification (GTV ) experiment 'was carried out 
in two distinct stages. The first stage involved modeling the ASTROMAST itself 
as it would be tested in the first open loop modal test, i.e., the beam alone in 
a cantilevered position. This produced modal frequencies, mode shapes, and mass 
integrals which were used in the second stage of the modeling process to deve- 
lop, through modal synthesis, a model of the entire GTV experiment including the 
AGS and the beam containment structure. . >• 

B. ASTROMAST Beam Model 

The ASTROMAST is a symmetric beam 512 inches in length and triangular of 
cross section. Three longerons form the corners of the beam and extpn-f along its 
full length unbroken. The cross members which give the beam its shape divide 
the beam into 91 sections having equal length and mass and similar elastic pro- 
perties. This sectioning was used advantageously in the modeling process. 

As a first step, a single section of the beam (See Figure 10.) was modeled 
from its geometry and the elastic properties of S-GLASS. This provided the 
necessary information for a 91 element lumped mass model which was used to 
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determine the modal frequencies, mode shapes, and mass integrals of the beam in 
a cantilevered position. Because of the syronetry of the beam, bending in the 
two transverse planes must be uncoupled and identical, i.e., bending modes in 
one transverse plane must be exactly repeated in the orthogonal transverse 
plane, and torsion nust also be uncoupled from these bending modes. This allowed 
a simplified analysis of the ASTROMAST beam via a model comprised of 91 infini- 
tesimally thin discs of mass connected by massless beam segments. A two dimen- 
sional model was used to find the modal frequencies, mode shapes, and mass 

integrals of the ASTROMAST beam in bending. In this model each of the 91 lumps 

0 

was allowed to translate in a direction transverse to the beam and to rotate 
about its own center of mass. (See Figure 11.) Also, shear and bending defor- 
mation were taken into account. Because the beam is very stiff longitudinally, 
deformation in this direction was considered negligible and, therefore, not 
modeled. The modal data resulting from the analysis of this model are shown in 
Table 1 along with the results of the torsional analysis. 

Again, because of the symmetry of the beam, the torsional analysis could be 
carried out separately from the bending analysis. The beam in torsion is a one 
dimensional problem because the lumped masses are allowed only to rotate about 
the longitudinal axis. If the beam is assumed to be uniform, a closed form solu- 
tion to the torsional bending problem can be found. This was done as a check of 
the numerical method; and indeed, the results of the two techniques compare 
favorably as shown in Table 2. 

C. Complete System Model 

Using the model of the ASTROMAST beam alone and a rigid body model of the 
AGS (Advanced Gimbal System), a complete system model was produced using modal 
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synthesis. The following are the major assumptions included in this model: 

1) the tip support structure has neglible effects upon the struc- 
tural behavior, 

2) the AGS is rigid, 

3) there is no stiffness associated with the translational motion 
of the base. 

The first assumption is easily justified by the fact that the tether 
"string'' which supports the tip will be approximately 60 feet long and very 
light as compared to the tip mass of about 32 pounds. Also, the tip deflections 
are expected to be at most a few inches in magnitude. The second assumption is 
justified when the stiffness of the AGS is compared to the relative flexibility 
of the A5TR0MAST beam. As for the third assumption, the base of the experiment 
would be unrestrained if it were not driven by the disturbance actuators. Of 
course the actuators will not only provide a disturbance input as desired but 
may also passively restrain the base in some way. This will be incorporated into 
the model as more data become available concerning the disturbance actuators. 
The mathematical model of the GTV experiment structure incorporates five rigid 
body degrees of freedom. Two of these result from the translations of the base; 
the other three from the rotations of the gimbal • system. In addition, the 
ASTROMAST is represented by its first seven flexible body modes to produce the 
system model summarized in Tables 3 and 4. 
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VI!. Modal Test of ASTROMAST Beam 


Modal tests have been conducted on the ASTROMAST beam at Marshall Space 
Flight Center. So that gravity would have as little effect as possible on the 
test results, the ASTROMAST was cantilevered in a hanging position for the modal 
test. In order to do this, a mounting bracket was fashioned as shown in Figure 
12 so that the AS T R0MAST could be mounted and then deployed downward unde r the 
influence of gravity. (See Figures 13 through 16.) 

When fully deployed, the ASTROMAST exhibits a longitudinal twist of about 
280 degrees as shown in the photograph of Figure 17. This twist contributes to 
coupling between the torsional and bending modes; and, although it was not 
includes in the model described in Section VI, the model seems to describe the 
beam quite well as is discussed later in this section. 

Testing using three different levels of input was found to produce the same 
experimental results with regard tc the modal model. This shows that the 
ASTROMAST behaves in a linear manner over a typical range of input excitations, 
thus making it an excellent choice of test articles for the LSS/GTV experiment. 

Typical modal damping for the ASTROMAST beam was found to be about 0.2 to 
0.3 percent. This is in the range perdicted for the ACOSS model #2 typical LSS 
which was constructed (on paper) of graphite epoxy. The ASTROMAST is similar in 
that the S-GLASS members from which it is constructed are also a composite of 
fiber and resin and would be expected to have similar damping characteristics. 

The lower modal frequencies of the cantilevered ASTROMAST beam as measured 
in the laboratory and as computed from the model described in Section VI are 
summarized in Table 5. The preliminary nature of the test data accounts for the 
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modes pedicted by the model which have not yet been identified and measured in 
the lab. However, it is clear that the model is in reasonable compliance with 
the actual structure. The modes determined experimentally thus far agree very 
well with the model as is revealed by the average percent error of about ten 
percent. Error in the model may be due to the unknown characteri sti cs of the 
S-GLASS as a composite. Since the actual percentage of S-GLASS present in the 
beam members was unknown, they were assumed to be made entirely of S-GLASS. This 
error in the modulus of elasticity could account for some of the modeling error 
present . 
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VIII. System Assembly, Verification, and Use 


The LSS/GTV experiment fixture will be assembled in the high bay area of 
building number 4419 at Marshall Space Flight Center. As mentioned earlier, each 
of the subsystems is being thoroughly tested independently in order to minimize 
the difficulties to be encountered upon assembly of the complete system. 

After assembly of the entire test apparatus, proper communication between 
the subsystems will be verified and proper operation of the subsystems in the 
test scenario will be established. The health status of each of the instrument 
packages will be of major concern at this time. Once each of the subsystems has 
proven to be operating properly in the test environment, the system will be 
ready to be tested in the operational configuration, i.e., to carry out a test 
of the first candidate control algorithm. 

Once the experiment fixture is shown to operate properly and consistently, . 
it will be available for testing on a regular basis and will be operated by 
Test Lab personnel. 
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IX. Summary and Future LSS/GTV Development 

MSFC is developing a LSS ground test facility to help verify LSS passive 
and active control theories. The facility is not parochial to the control area 
alone but has a wide purview of activities in which experimentation can be 
accomplished. These areas include but should not be limited to the following: 

1. Subsystem and component testing, 

2. Remote sensing and control, 

3. Parameter estimation and model verification and, 

4. Evolutionary modeling and control. 

All of these areas are currently in the planning stages of the facility develop- 
ment and should be an integral part of the facility within the immediate 
future. 
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MODE 

FREQUENCY (HZ) 

FIRST BENDING MODE 

0.613 

SECONO BENDING MODE 

3.917 

THIRD BENDING MODE 

11.179 

FIRST TORSIONAL MODE 

6.877 

SECOND TORSIONAL MODE 

20.628 

THIRO TORSIONAL MODE 

34.373 

■ - J 


Table 1. Model Frequencies of Cantilevered ASTROMAST Beam 
Model in Bending and Torsion. 


r ~ ~ ~ ~ - . ^ 

MODE 

CLOSED FORM FREQUENCY (HZ) 

TORSIONAL BEAM FREQUENCY (HZ) 

FIRST TORSIONAL MODE 

6.92 

6.83 

SECOND TORSIONAL MODE 

20.74 

20.63 

THIRD TORSIONAL MODE 

34.53 

34.37 

L . . 


Table 2. Torsional Mode Frequencies Resulting from Analytical Uniform Beam 
Closed Form Solution and 91 Element Lumped Mass Numerical Solution. 
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Table 4. LSS/CTV Experiment Structure Mathematical Model Flexible Body Modes. 





* Two modes at each frequency 

I 

table 5. Summary of Modal Frequencies from Model and Test Results for the 
Cantilivered ASTROMAST Beam. 
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APPENDIX H 


ATTACHMENT FROM MR # 23 REVIEW OF TECHNICAL PAPERS 



ATTACHMENT A. Review of Technical Papers 


A.l Introduction 

A total of twelve (12) memos and papers by Dr. Henry Waites were supplied to 
Control Dynamics by NASA for the company's information and use in the study of 
large space structures. Ten (10) of these papers deal directly with control 
systems, such as observer designs or closed-loop pole placement methods. Of the 
ten papers there are three principle concepts treated: observers, closed- loop 

pole placement and a disturbance isolation technique. 

Three of the ten papers were selected for critical review as they embodied 
the three basic concepts. The objectives of the review were: 

(1) Check and verify the equations and derivations, 

(2) Relate these new techniques to "standards” in the literature, 

(3) Identify strengths and weaknesses of the methods, 

(4) Determine suitable topics for further study using these methods. 

A. 2 Observer Techniques 

Two Ovservers were compared in the first section: Dr. Waites' Technique as 

described in the AIAA paper "An Observer for a Deployable Antenna. "[1] and the 
classic Luenberger[2] observer. Their block diagrams are depicted in Figures 1 
and 2, respectively. 

Dr. Waites' model circumvents the problem usually encountered with obser- 
vers, the inputs (particularly the disturbances must constitute part of the 
input to the observer as can be seen in Figure 2. The Waites' model proposes to 
utilize the rate-of- the- state signals (as measured by accelerometers) that 
include the vehicle responses to all disturbances as well as the command inputs. 
Although this idea does indeed eliminate the need to know the characteristics of 
the disturbance inputs, it also imposes an additional constraint matrix which 
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may result in a more sensitive system. The main advantage of the Waites' 
observer is that it provides a means of incorporating the disturbance inputs 
into the observer by direct measurement. 

If time is available for future work, a topic for additional study is the 
sensitivity of the observer system to variations in sensor locations from the 
point of disturbance inputs. 

A. 3 Closed-Loop Pole Placement Methods 

This section discusses two methods developed by Dr. Waites in memorandum 
ED/2-81-1, "Control Pole Placement Method." One is designated a control pole 
placement method with specified filters and the second is called a 
filter/control pole placement technique. The equations used by Waites can be 
expressed in block diagram form as shown in Figure 3. This one diagram is 
applicable to both methods. In the first method only the gain matrix F is 

unknown and is determined by Waites' method. The second technique requires the 

determination of both the gain matrix F and the filter matrix D. 

The first step in both methods Is to solve the closed-loop system equation 
of Figure 3 for the characteristic equation, given in matrix form as 

_ l 

AU) = | Xl m - A 1 1 Xl p - D - EC(\I m - A) BF | (10 a ) 

In method I matrices A, B, and C are known and matrices D and E are speci- 

fied; that is the designer has to select the elements of these two matrices. 
The elements of D should be distinct from the values of the A matrix. Method II 
starts with the equation above but only the E matrix is specified. The second 
determinant of ( 10 a ) 

_i ’ 

A (X i ) = | Xi I p - D - EC ( X i I m - A) BF | 
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Figure 3. WAITES POLE PLACEMENT METHODS 

SYSTEM EQUATIONS 
x = Ax + Bu 
y = Cx 

4 = Da + Ey • 

<jl= Fa 
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is solved for any row, by equating that row to zero. This will require the de- 
signer to select a maximum of r + p Ai's in order to solve for the resulting 
elements of D and F. As with almost any design technique, in order to use these 
methods a considerable amount of experience with the methods will be required. 
Also a solution for matrices D and F does not necessarily yield a state closed- 
loop system, as can be seen from Waites' simple example in the paper which has 
the unspecified closed-loop pole in the right half-plane. 

Dr. Waites' derivation was examined extensively and is correct. As can be 
seen from the block diagram of Figure 3, the system equations yield a typical 
system with a filter and a couple of gain matrices in the feedback loop. In 
method I after the selection of a specific number of closed-loop poles the open- 
loop poles of the filter are selected arbitrarily and then solved for the gain 
matrix F. 

In method II the closed-loop poles are selected and then the open-loop 
poles of the filter are found along with the elements of the gain matrix F. 
Some questions remain regarding this method, such as gain values, performance, 
and stability characteristics. These methods may provide techniques that can be 
used to rapidly establish a satisfactory control system since they were devel- 
oped specifically for computer implementation; but as with most new design tech- 
niques, more experience using the methods may be necessary. An obvious area for 
further study of these techniques would be a computer oriented procedure to map 
the possible boundaries of the unassigned closed-loop poles as a function of the 
assigned poles. With such a procedure, more information would then be available 
regarding system stability and performance. 

A. 4 Disturbance Isolation Controller 

The controller described in memorandum ED/12-81-5, "Disturbance Isolation 
Control 1 er" , is another application of observer theory. In developing the 


H-5 


example given in the memo. Dr. Waites defines a generic model. This model is 
shown with solid lines in Figure 4. The observer of Figure 4 appears to be a 
combination of a Luenberger and modified Waites observers. The dashed lines in 
the figure represent one possible implementation of the disturbance isolation 
feature for the example model. 

An example is given in the memo to illustrate the control scheme. This 
example starts out with an eighteenth order unobservable system. By manipu- 
lating the equations the system is reduced to a twelve order observable system 
and arranged in state variable format. 

To implement the disturbance isolation feature, three gain matrices are 
placed in the feedback loops as shown with dashed lines in Figure 4. These gain 
matrices are shown as K, U and T, where: 

K = [0 0 -Ko -Kj] 

U = [0 (m32 - 11131 dsm) 0 0], and 

Pj = C m 31 "SI dsm 0 ] 

(these symbols K, U and Pj are not used in Waites memo). 

Waites uses the reduced order example model to develop the observer then he 
returns to the full -order model to demonstrate the disturbance isolation feature 
and this raises a question regarding the disturbance vector F. As can be seen 
from Figure 4, the input Z to the observer is the sum of Px and Pq F. The x 
signals are assumed to be measured with accelerometers but there is no apparent 
assumed method for measuring F. Regular disturbances like thruster firings can 
be characterized, but random disturbances such as astronaut movements are more 
difficult to model. 

An area for possible further study would be to define a suitable system 
model, develop a disturbance isolation controller, and then conduct a system 
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Figure 4 DISTURBANCE ISOLATION CONTROLLER 
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performance evaluation by characterizing the disturbance function with known and 
unknown components. 

A. 5 Summary 

A characteristic of the methods reviewed that have not been discussed is 
the robustness of the techniques. That is, how will the closed-loop system per- 
form for the expected uncertainties in the parameters of the system's state 
equations? The only way to answer such a question is with realistic case stu- 
dies. Because of the apparent complexities of the Waites observer and the 
disturbance isolation controller, it is expected that these two methods will not 
exhibit as much robustness as the pole placement method. 
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GROUND TEST experiment for urge space structures 
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ABSTRACT 


In recent years a new body of control theory 
has been developed for the design of control 
systems for Large Space Structures (LSS). The 
problems of testing this theory on LSS hardware are 
aggravated by the expense and risk of actual "in 
orbit" tests. Ground tests on large space 
structures can provide a proving ground for 
candidate control systems, but such tests require a 
unique facility for their execution. The current 
development of such a facility at the NASA Marshall 
Space Flight Center (MSFC) is the subject of this 
paper. 


INTRODUCTION 

As the U.S. moves Into the Shuttle era of 
space technology, there are numerous proposals from 
the scientific, civilian ano Defense Communities 
which envision the use of Large Space Structures 
(LSS). By definition, a LSS Is very flexible, 
lightly damped, and exhibits multiple vibrational 
modes of very low frequency. Many of the missions 
alluded to above require high performance from the 
LSS in areas such as precision pointing and 
preservation of vibration free Image planes. To 
test the new and various schemes proposed for LSS 
control, Marshall Space Flight Center (MSFC) Is 
establishing a LSS laboratory In which 
experimentation with LSS-like structures may be 
performed. The first experiment to be Installed in 
the laboratory Is being developed concurrently. 


SYSTEM OVERVIEW 

The Ground Test Verification (GTV) experiment 
Is described by the drawing of Figure 1. The first 
test article will be the ASTROMAST beam as shown. 
The ASTROMAST is extremely lightweight (about 5 
pounds) and approximately 45 feet In length and is 
constructed almost entirely of S-GLASS. 

The test article will be mounted on the 
faceplate of the Advanced Gimbal System (AGS) 
engineering model Uiich, along with an additional 
torque actuator In the azimuth, provides the 
control inputs for the system. The azimuth gimbal 
also provides a means of rotating the entire 
experiment manually to produce different test 
scenarios. The ASTROMAST will be gravity unloaded 
by a cord extending downward from the base 
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excitation system of the ASTROMAST to the tip. The 
tension in this cord will be adjustable to provide 
a means of "tuning" the structure. 

The AGS will be supported by the base 
excitation system of the beam containment structure 
which will be free to translate In the horizontal 
plane and will Include hydraulic actuators to 
provide translational disturbance Inputs to the 
test fixture. These disturbances will represent 
Astronaut push-off or RCS (Reaction Control System) 
thruster firing. 

Six separately packaged Inertial measurement 
assemblies comprise the control system sensors. Two 
of the packages, containing three axis 
translational accelerometers, are Identical. One 
will be mounted on the mast tip, and the other on 
the AGS side of the base excitation system. Three 
other packages contain ATM (Apollo Telescope Mount) 
rate gyros and will be Installed on the AGS 
faceplate. The sixth package, the Kearfott Attitude 
Reference System (KARS), will be placed at the mast 
tip along with the accelerometer package. 

The signals from these instruments will be 
read by the COSMcC I data gathering and control 
system, and processed according to the control 
strategy under scrutiny. The control actuator 
signals will then be transmitted to the AGS as 
Inputs to the «iynamical system. 

The COSMEC I will be Interfaced to a Hewlett 
Packard HP9845C desktop computer which will store 
data as It Is collected during a test run, and then 
provide post experiment data reduction and display 
off-line. The controller inputs and outputs 
(measurements and commands) can be recorded at each 
sample period or at some multiple of sample 
periods. 


SUBSYSTEMS 

The subsystems which comprise the LSS/GTV 
experiment fixture as described are currently in 
various stages of development at NASA MSFC. 
Thorough verification of each of the subsystems in 
controlled test environments comprises a 
significant part of the preliminary system testing. 
The subsystems will not only be tested 
individually, but will be tested in an Integrated 
laboratory environment where each of 
the subsystems will interact with the others in a 


1-1 



wanner much If an actual test article were being 
used. Such testing Is designed t ensure proper 
operation of the complete test fixture upon 
assembly. 

ADVANCED GIMBAL SYSTEM 

The Advanced Gtmbal System (ACS) Is a 
precision, two axis glmbal system designed for high 
accuracy pointing applications. The AGS gimbals 
serve the elevation plane and a third glmbal has 
been added to the system In the azimuth. The AGS 
receives torque commands from the COSMEC I data and 
control system In the form of analog Inputs over 
the range of -10 to +10 volts. Because the AGS 
servo amplifier outputs a current which causes an 
applied torque proportional to the current, the 
control algorithms used In the COSMEC 1 must be 
designed to produce torque command signals. The 
AGS glmbal torquers can generate 37.5 ft. -lbs. of 
torque and the azimuth torquer can generate 13.8 
ft. -lbs. 

COMPUTATIONAL ELECTRONICS 

The COSMEC I Is an AIM 65 based (MC6502 
processor) micro system which is used to handle 
data from the control system sensors, output 
commands to the control system actuators, transmit 
data for storage to the HP 9845C desk-top computer, 
and Implement the control and Inertial strapdown 
algorithms. The COSMEC I uses special hardware and 
software to allow the handling of a variety of 
devices (sensors, actuators, etc.) in real time. It 
also makes use of four hardware arithmetic 
processors to reduce computation time. 

The COS.MEC I "reads" a variety of types of 
sensor output signals via Interface cards which are 
an Integral part of the COSMEC I system. These 
cards allow the COSMEC I processor to Interface 
with each of Its peripherals In a similar manner. 

The software used In the COSMEC I system may 
be separated Into four basic groups: 1) utility 
software for handling the various hardware cards 
which Interface to Instruments, 2) software to 
Implement the control algorithm, 3) software to 
Implement the Inertial strapdown algorithm, and 4) 
Initialization and startup software to ready the 
Instruments and equipment for a test. 

The digital stroller software for the first 
ground test experi.r*.it will Implement a linear 
discrete multivariable controller having multiple 
Inputs and outputs. The controller will be In state 
variable form and will be programme so that the 
system matrices are Initial Input data to the 
program and can be stored on tape and easily 
changed. 

Because the inertial measurement Instruments 
measure with respect to Inertial reference space, 
there Is a natural bias In the measurements due to 
the acceleration of gravltv and earth's rotation. 
That Is, In the earth based experiment the 
accelerometers measure about one g acceleration 
downward and the rate gyros measure about 15 deg/hr 
rotation wnile at rest with respect to the 


laboratory reference frame. The Inertial strapdown 
algorithm provides a means of removing this bias 
from the miasurement Instruments. 


INERTIAL MEASUREMENT ASSEMBLIES 

Three different types of Inertial measurement 
assemblies are planned for use on the Ground Test 
Verification structure In the first experiment: the 
Kearfott Attitude Reference System (KARS), the 
Apollo Telescope Mount (ATM) rate gyros, and two 
accelerometer packages developed by NASA. Each of 
the Instrument packages generates signals In a 
particular form different from the other 
Instruments as was mentioned In the section dealing 
with the COSMEC I Interface cards. These different 
signal types are discussed In the following as each 
of the Instruments Is discussed Individually. 

The Kearfott Attitude Reference System (KARS) 
Is an attitude measurement system designed for use 
In the U.S. Army remotely piloted vehicle. It 
provides measurement resolution of 13.9 X 10*3 
deg/sec In the pitch and yaw axes (axes transverse 
to the ASTROMAST) and 25.0 X 10*3 deg/sec In the 
roll axis (ax<s along the length of the ASTROMAST). 
The dynamic range of the rate gyro outputs of the 
KARS Is 40 deg/sec In pitch and yaw and 70 deg/sec 
In roll. Because of Its light weight, (8.9 pounds) 
the KARS will be used as the mast tip rotation 
sensor In the first ground test experiment. 

The output signals of the KARS are in the form 
of asynchronous digital pulses. One signal, the 
change In angular position In yaw for Instance, 
requires two channels; one for pulses representing 
positive rotation and .the other for pulses 
representing negative rotation. The COSMEC I 
system accumulates the pulses over a 20 millisecond 
period to produce measurements of the angular rate 
and position of the ASTROMAST tip. 

The Apollo Telescope Mount (ATM) rate gyro 
packages are designed to measure small angular 
rates very precisely. Each package measures angular 
rate In one axis with resolution finer than 0.5 X 
10* 3 deg/sec and offers J dynamic range of i 1.0 
deg/set. The ATM :te gyro packages will be mounted 
on the faceplate of the engineering AGS so that 
they will measure the rotation of the base of the 
test article. 

The output signals of the ATM rate gyro 
packages are 145 volt analogs and ere handled by 
the analog to digital converter card of the COSMEC 
I system where they are converted to 12 bit binary 
words. 

Two Identical accelerometer packages will be 
used on the ground test experiment fixture in the 
first test. One package will be placed on the mast 
tip along with the KARS and the other on the test 
fixture base excitation system as shown In Figure 
1. The accelerometers provide resolution finer 
than .0001 g and a :ynanic range of t 3 g with a 
bandwidth of 25 to 30 Hz. 

The signals from the accelerometers are 
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different from either those of the KARS or the ATM 
rate gyros. As In the cue of the KARS, two 
channels ere required for each of the degrees of 
freedom of the iccelerometer package, l.e., si* 
channels per accelerometer package. One channel of 
each pair carries a 2.4 kHz. square wave 
Synchronization signal and the other channel 
carries the acceleration information. Zero 
acceleration is represented by a signal Identical 
to that of the synchronization channel , positive 
acceleration by an Increase In frequency, and 
negative acceleration by a decrease In frequency as 
compared to the synchronization channel. As In the 
rases of the other Instruments, these signals are 
monitored by a hardware card In the COSMEC I 
system. 

BEAM CONTAINMENT STRUCTURE 

The Oeam containment structure Includes the 
base excitation system, the disturbance actuators 
and signal source(s), and power supply for the 
entire test fixture. This Is essentially that 
equipment required of a laboratory to carry out 

dynamic tes ting of structures such as the ground 

test experiment. As depicted In Figure 1, the beam 
containment area can accommodate structures 
approaching 120 feet In height. Access Is provided 
at various levels along the structure via catwalks. 
Also the control room for experimental operations 
Is at a position about 50 feet above the floor of 
the test facility, thus making the experiment 
easily accessible. 

The base excitation system for the GTV 
experiment Incorporates a linear bearing 
arrangement which is designed to restrain the 
experiment as little as possible In the two 
translational degrees of freedo.a while allowing no 
rotation. Translation of the base of the AGS via 
the base excitation sy'-'m provides the means of 
applying disturbances the system. The base 
excitation system uses ilgh pressure hydraulic 
system to effect force ...,.uts to the structure. 


SYSTEM MODELING AND PERFORMANCE 

For purposes of system studies and controller 
design, an analytical model of the Ground Test 
Experiment Is necessary. Modeling of the Ground 
Test Verification (GTV) experiment was carried out 
In two distinct stages. The first stage involved 
modeling the ASTRCMAST Itself a: It would be tested 
In the first open loop modal test, l.e., the beam 
alone In a cantilevered position. (See Table 1 for 
a listing of the resulting modal frequencies.) This 
produced modal frequencies, mode shapes, and mass 
Integrals which were used In the second stage of 
the modeling process to develop, through modal 
synthesis, a model of the entire GTV experiment 
Including the AGS and the beam containment 
structure. The modal frequencies resulting from 
the analytical model of the complete experiment 
appear In Table 2. 

To date, modal tests have been conducted on 
the ASTRCMAST bean alone In a cantilevered 


conf tguratlon at Marshall Space Flight Center. So 
that gravity would have as little effect as 
possible on the test results, the ASTROMAST was 
cantilevered In a hanging position. When fully 
deployed, the ASTROMAST exhibits a longitudinal 
twist of about 230 degrees which contributes to 
coupling between the torsional and bending modes. 

Measured mndal frequencies resulting from the 
tests are shown In Table 3 along with the predicted 
mudal frequencies and percentage errors. The 
errors are reasonable given the limited amount of 
knowledge about the percent composition of the S- 
GIASS composite and the unexpected twist which was 
not Included In the mathematical model. 


SUMMARY AMO FUTURE DEVELOPMENT 

Marshall space Flight Center Is developing a 
Large Space Structure Ground Test Verification 
experiment facility having adequate fidelity and 
flexibility to accommodate the demands of LSS 
control theory testing. The first experiment is in 
the subsystem verification and Integration phase 
with the first "all up" test targeted for spring 
1984. This test employs the ASTRCMAST, a 
lightweight S-GLASS composite deployable beam 
structure, as the test article and is cited to 
prove out centralized and distributed sensor 
control strategies. 

Future plans for the facility Include the test 
of more complicated structures and sensor/actuator 
arrangements. The first of these will be the 
addition of structural hardware to the existing 
ASTROMAST In order to produce a rcructure having 
more LSS- T 1 ke characteristics. 


MODE 

FREQUENCY (HZ) 

FIRST BENDING MODE 

0.618 

SECONO BE NO INC MODE 

3.91? 

THIRD BENOING MOOE 

11.179 

FIRST TORSIONAL MOOE 

6.877 

SECONO TORSIONAL MOOE 

20.628 

THIRO TORSIONAL MOOE 

34.373 


Table 1. Modal Frequencies of Cantilevered 
ASTRCMAST Beam Modei In Bending 
and Torsion. 
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